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ABSTRACT 


AN  EXPERIMENTAL  STUDY  OF  THE  TRANSPORT  OF  A 
NON-DIFFUSIVE  SCALAR  CONTAMINANT 
IN  THE  DECAYING  TURBULENCE  FIELD 
OF  AN  ENCLOSED  CHAMBER 


An  experimental  investigation  of  the  convective  properties 
of  a turbulence  field  within  a closed  mixing  chamber  has  been 
executed;  the  investigation  is  supported  by  appropriate  theo- 
retical analyses.  The  purposes  of  the  study  were  (i)  to  develop 
an  experimental  technique  which  could  ascertain  the  convective, 
or  dispersive,  properties  of  such  a turbulence  field  (i.e.,  with 
the  constraint  of  "zero"  molecular  diffusivity)  and  (ii)  to 
document  the  influence  of  the  Reynolds  number  on  the  convective 
transport  properties  of  a representative  field. 

The  derived  technique  involves  the  introduction  of  sub- 
micron particulate  matter  into  the  lower  half-volume  (as  defined 
by  a horizontal  splitter  plate  which  is  rapidly  removed  from 
the  chamber  at  the  initiation  of  an  experiment)  and  the  sub- 
sequent monitoring  of  the  presence/absence  of  this  contaminant 
along  a vertical  line  centered  at  the  plane  of  the  splitter 
plate.  An  optical/mechanical/electrical  technique  was  used  to 
scan  205  locations  along  a 7.4  cm  portion  of  an  incident,  focused 
(d  = 0.1  mm)  , laser  beam  and  the  scattered  light  intensity  was 
used  to  identify  the  presence  or  absence  of  the  fluid  elements 
which  originated  in  the  lower  half-volume.  That  is,  since  the 
particulate  matter  is  shown  to  be  too  large  to  respond  to  the 
diffusion  effects  of  the  Browman  motion  of  the  gas,  and  too 
small  to  move  independently  of  the  fluid  element  of  its  origin, 
it  is  inferred  that  it  faithfully  marks  the  lower  half-volume 

xi  -J9-  
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fluid  over  the  duration  of  the  experiment  (0.6  sec).  The  7.4  cm 

scan  length  covered  2.92  mesh  lengths  of  the  biplane,  round  rod, 

turbulence  generating  grid.  Air  was  used  for  69  experiments; 

Freon  12  was  used  for  100  experiments  (v  . = 5.8v„  , _). 

r air  Freon  12 

A principal  result  of  this  investigation  is  that  the 
particle  dispersion  is  dependent  upon  the  mesh  Reynolds  number 
(UqM/v) ; the  width  of  the  mixing  field  was  found  to  decrease 
as  the  Reynolds  number  was  increased.  A theoretical  model  of 
the  dispersion  was  developed  using  the  suggestion  by  Batchelor 
and  Townsend  [1956]  that  the  Lagrangian  velocity  auto  corre- 
lation function  can  be  made  stationary,  in  a decaying  turbulence 
field, by  an  appropriate  stretching  of  the  time  coordinate.  TK? 
present  model  utilizes  both  the  Lagrangian  micro  and  macro 
scales  to  describe  this  correlation  function;  numerical  values 
for  these  two  scales  are  inferred  from  the  experimental  data. 

The  observed  decrease  in  the  particle  dispersion  with  increasing 
Reynolds  number  is  supported  by  a previously  suggested  re- 
lationship between  the  Lagrangian  micro  scale  and  the  mesh 
Reynolds  number,  a new  value  for  the  coefficient  in  this  re- 
lationship is  proposed. 

Analytical  considerations  are  presented  to  support  the  use 
of  a (spatial)  Gaussian  cumulative  distribution  function  for 
the  description  of  the  ensemble  average  concentration;  viz, 

T (z,t)  = [ 2 tt  am  ( t ) 2 ] ~ 1/2  I2  exp  _ [(fi  - <z>)/2om(t)2]dd. 
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Satisfactory  agreement  between  this  distribution  and  the  experi- 

2 

mental  data  was  obtained;  ~ 


°m(t> 

square  particle  dispersion,  z^ 


was  used  to  evaluate  the  mean 
Two  properties  of  the  instan- 
taneous scans  were  defined:  i)  an  ensemble  average  measure  of 
the  instantaneous  width  of  the  mixing  field,  o^,  and  an  ensemble 
average  measure  of  the  migration  of  the  instantaneous  center 


lines,  o 


Az ' 


The  instantaneous  scans  are  shown  to  have  essentially 


the  same  width  as  that  of  the  ensemble  average  field  o. 


0 . 8a 


whereas  exhibits  very  little  growth  with  respect  to  time 

following  the  establishment  of  its  nominal  value  at  U t/M  - 8. 

o 


m 


xii 


NOMENCLATURE 


X 


B 

B (k) 


B* 


radius  of  particle 
a.  gstrom  10  ^ 

area  of  exit  discharge  nozzle 
experimental  constant  in  energy  decay  law, 
area  of  piston 

cross  sectional  area  of  grid 
constant 

Fourier  transform  of  !'(x,y,z) 
complex  conjugate  of  B 

2 

amplitude  spectra  of  < •»  > 

concentration  of  scalar  contaminant 

constant  7.14  cm/radian 

constant  defined  in  (6.3.1) 

drag  coefficient 

strength  of  irrotational  vortex 

see  (3.2.24) 

Cx/(!-n/2) 
see  3.2.5 
yil-n/2) 

diameter  of  Diplanar  grid  elements 

minimum  diameter  of  focused  laser  beam 

1 2 

diameter  of  laser  beam  at  ( /e  ) intensity 
molecular  diffusivity 
turbulent  diffusivity 
particle  diffusivity 
focal  length  of  lens 


see  3 . 
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drag  force  of  grid 
Fourier  transform  of  u(x,y,z) 
gravitational  constant 
lateral  dimension  of  chamber 

distance  between  concentrated  vorticity  centers 

intensity  of  segregation 

indicial  notation 

Boltzman  constant 

wavenumber  variable 

length  of  cylinder  (hydraulic) 

scale  of  segregation 

length  scale  for  the  dispersion  ar^lysis,  see  (3.2.21) 
mass  of  particle 

total  mass  which  is  attached  to  piston  rod 
mesh  length  of  grid 

m = (l-n/2),  m is  used  in  the  dispersion  analysis 

number  of  experiments 

exponent  of  decay  law,  see  (3.1.7) 

pressure  (static) 

pressure  acting  on  piston 


-2  -2  -2 

turbulent  energy  intensity  (u  +v  +w  ) 


1/2 


radial  coordinate 

inner  radius  of  helical  slit  6.99  cm 
eddy  of  Kolmogoroff  size 

Lagrangian  correlation  coefficient,  see  (3.2.19) 
Reynolds  number  (UoM/v) 

Schmidt  Number  (v/D) 

Strouhal  number 


xiv 


time  variable 

time  increment  equal  to  the  transform,  (3.2.24),  of 

the  microscale  tt 
L 

apparent  origin  of  kinetic  energy  decay 


Tn(k) 


from  the  transformation  of  the  juncture  point,  see  (3.2.25) 
temperature  (absolute) 

total  elapsed  time  for  passage  of  grid  through  chamber 
see  (3.2.21),  release  time  of  the  contaminant 
transfer  function 
velocity  in  x direction 

component  of  velocity  (cylindrical  coordinates) 

speed  of  grid 

velocity  in  y direction 


radial  component  of  velocity  (cylindrical  coordinates) 
velocity  in  z direction 


<Xi2> 


r.m.s.  of  the  z component  velocity  fluctuation 
total  work  done  on  fluid  by  action  of  grid 
mean  square  displacement  of  a fluid  particle 


(x,y,z)  cartesian  coordinates,  x in  direction  of  U , z in 
direction  of  mean  gradient  of  r 


2 

<z  > 

P 


length  of  focused  laser  beam  in  which  d is  constant 

o 

centerline  of  mixing  chamber 

normalized  root  mean  square  of  the  particle  displacement 

z 2/L2  = z*2 
p p 

center  of  ensemble  average  distribution 

mean  square  displacement  of  particle  in  the  z direction 
see  figure  29  for  definition  of  special  symbols 
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dummy  variable  of  integration 
3p/(2pp+p) 

concentration  of  non  diffusive  contaminant 

fluctuating  value  of  I';  ( r — < r > ) 

incremental  value 

boundary  layer  thickness 

dissipation  of  turbulent  energy 

frequency  (wavenumber ) 

radians  (angular) 

unit  vector  in  0 direction 

wavelength  of  disturbance 

Taylor  microscale  (g:  spatial,  t:  time) 

wavelength  of  laser  radiation 
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absolute  viscosity 
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density  of  fluid 

density  of  particle 
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<r  (z,t)  > (a  Gaussian  distribution  for  ^l'(7.,t)s  is 
assumed,  see  section  5)  , see  3.2.12 

see  5.4.1 

see  5.4.2 

stretched  time  coordinate,  i * C In  At/At^,  see  (3. 2.. ’4, 5 
Lagrangian  microscale  in  the  stationary  coordinates 
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Lagrangian  velocity 
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1.  INTRODUCTION 

The  efficient  mixinq  of  one  fluid,  or  a contaminant  sub- 
stance, with  the  host  fluid  is  an  intrinsic  feature  of  a tur- 
bulent flow.  The  term  "mixinq"  is  often  used  to  describe  the 
net  effects  of  turbulent  dispersion,  in  which  the  eddyrnq 
motions  provide  a convective  transport  of  the  material  to  be 
mixed,  and  molecular  diffusion,  in  which  the  concentration  of 
the  material  in  a given  fluid  element  is  changed  by  the  trans- 
port associated  with  the  aqitation  at  the  molecular  level.  As 
noted  by  Brodkey  [1975)  the  interpretation  of  the  deqree  of 
mixing  is  related  to  the  scale  of  the  observation,  i.e.,  to 
the  size  of  the  volume  used  to  define  the  concentration.  Two 
measures  of  mixing  are  commonly  employed,  Brodkey  [1975],  the 
intensity  of  segregation  I and  the  scale  of  segregation  Lg . 

For  a binary  mixture,  where  C is  the  instantaneous  concentration 
(C+c)  of  the  contaminant,  the  intensity  of  the  segregation 
(1^)  is  given  as 

2 

1 = ~ 

S C(l-C)  (1.1) 

and  the  scale  of  the  segregation  (I,  ) is  given  as 

L = / c(x)c(x+r)dr  (1.2) 

The  turbulent  motions  play  a dual  role  in  the  mixing  pro- 
cess. They  are  responsible  for  the  convective  transport  of  the 
material  to  be  mixed  and  they  can  greatly  enhance  the  molecular 
diffusion  by  straining  the  concentration  field  and  hence  steepening 
the  spatial  gradient  of  the  transported  material.  This  dual  role, 
and  the  influence  of  the  initial  events  on  the  subsequent  dis- 
tributions, result  in  turbulent  mixing  phenomena  which  are  of 
fundamental  interest. 

The  present  investigation  was  conceived  and  executed  as 
the  first  stage  of  a research  program  to  examine  these  basic 
phenomena  of  turbulent  mixing.  The  total  program  is  to  involve 
measurements  of  the  velocity,  the  mol  ^ular  species  concentration 


1 


2 


and  the  presence  or  absence  of  a non-di f f us ive  contaminant  in 
a simple  mixing  field  which  may  be  characterized  as  a decaying, 
homogeneous  turbulence  field.  The  initial  condition  for  the 
mixing  is  similarly  simple;  the  fluids  to  be  mixed  are  segregated 
in  two  distinct  volumes  until  the  initiation  of  the  mixing 
process.  The  series  of  experiments  in  the  total  program  are 
to  reveal  the  influence  of  the  control  variables:  different 
states  of  the  initial  turbulence  kinetic  energy  and  scales, 
and  ii)  different  densities  of  the  gases  to  bo  mixed.  The  data 
acquisition  techniques,  discussed  below,  are  to  provide  quite 
detailed  representations  of  the  mixing  process. 

The  present  study  is  solely  concerned  with  the  developmei  i 
and  application  of  the  non-di f f us ive  contaminant  experimental 
technique.  The  contaminant  is  introduced  into  one-half  of  the 
mixing  chamber  and  is  convectively  mixed  with  the  other  half- 
volume following  the  initiation  of  the  experiment.  Mechanical 
optical  techniques  are  used  to  identify  the  presence/absence  oi 
the  contaminant  within  small  volumes  along  a segment  of  the 
* mixing  region.  Since  the  half-volume  of  origin  of  the  fluid 

element  is  known,  the  dispersive  properties  of  the  turbulence 
field  may  be  inferred.  The  use  of  a closed  mixing  chamber 
greatly  simplifies  the  task  of  providing  a sufficient  volume 
of  the  marking  contaminant.  Hence,  this  distinctive  experimen 
technique  would  be  somewhat  more  difficult  to  utilize  in  con- 
ventional wind  tunnel  flow  systems. 

The  first  analytical  description  of  dispersion  in  a 
homogeneous,  stationary  turbulence  field  is  that  provided  by 
Taylor  (1921)  in  the  classic  paper  entitled  "Diffusion  by 
Continuous  Movements,"  Townsend  (19r)4)  recognized  that  the 
straining  action  of  the  turbulence  would  influence  the  (molecular) 
diffusive  as  well  as  the  dispersive  properties  of  the  field 
and  suggested  appropriate  modifications  to  the  Taylor  analysis. 
Townsend's  concepts  wore  refined  by  Saffman  (1961)  and  the 
full  problem  has  been  reconsidered  most  recently  by  Chevray 
(1977) . 
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The  physically  realistic  case  of  decaying  homogeneous 
turbulence  is,  in  contradistinction  to  the  above  studies, 
non-stationary . Batchelor  and  Townsend  (1956)  propose  that  an 
appropriate  stretching  of  the  time  coordinate  will  cause  the 
Lagrangian  correlation  function  to  be  stationary.  This  hypothesis 
is  a central  element  of  the  theoretical  structure  which  is 
utilized  to  interpret  the  experimental  results.  The  objective 
of  these  theoretical  considerations  is,  like  the  prior  efforts 
of  Uberoi  and  Corrsin  (1953)  and  Shlien  and  Corrsin  (1975),  to 
deduce  Lagrangian  scales  from  the  observed  dispersion  measure- 
ments. Tennekes  [1975]  has  reexamined  the  data  of  Shlien  and 
Corrsin  and  has  derived  an  analytical  relationship  between  the 
Eulerian  and  Lagrangian  time  microscales  in  decaying  isotropic 
turbulence.  Tennekes  shows  that  the  ratio  of  the  Lagrangian 
microscale  to  the  Eulerian  time  microscale  is  proportional  to 
R^;  ( R ^ = q A/v)  . Quantitative  reevaluations  of  the  Shlien 
and  Corrsin  and  the  Tennekes  results  are  presented  in  Appendix  E. 
An  important  difference  between  the  present  and  former 
experiments  is  the  diffusive  nature  of  the  contaminant  used  to 
mark  the  fluid;  a line  source  of  heat  was  used  to  thermally  mark 
the  fluid  elements  in  the  above  two  research  efforts.  A second 
major  difference  is  associated  with  the  time  which  is  available 
for  the  dispersion.  The  temperature  "tagging"  experiments 
were  initiated  following  the  establishment  of  a homogeneous 
and  approximately  isotropic  turbulence  state.  Consequently,  a 
relatively  limited  extent  of  the  stretched  time  coordinate  was 
available  for  the  contaminant  dispersion.  Conversely,  the 
present  experiment  initiates  the  dispersion  in  the  initially 
inhomogeneous  field  and  thereby  gains  a substantial  period  of 
development  in  the  stretched  time  coordinate.  A third  feature 
of  the  present  experiment  is  that  the  primary  data  are  derived 
from  essentially  instantaneous  scans  through  the  mixing  field 
with  the  ensemble  statistics  being  generated  from  a number 
of  experiments.  Consequently,  the  data  base  is  characterized 
by  considerable  detail  in  terms  of  a given  scan  through  the 
mixing  region  and  a relatively  limited  number  of  realizations 
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with  which  to  form  any  particular  ensemble  average.  The  former 
feature  is  used  to  evaluate  the  statistical  features  of  the 
instantaneous  scans;  these  interpretations  are  not  possible 
with  the  classical  wind-tunnel,  single-probe  techniques. 

The  above  matters  are  made  explicit  in  the  following 
sections.  Section  2 presents  the  Experimental  Facility  and 
Techniques,  including  the  mixing  chamber,  the  optical  data 
acquisition  system,  the  contaminant  generator  and  the  signal 
processing  units.  Analytical  descriptions  of  the  various  aspects 
of  this  investigation  are  presented  in  Section  3.  The  dynamics 
of  small  particles  and  the  analytical  considerations  which 
support  the  interpretation  of  the  present  particles  as  being 
non-dif fusive  are  presented  in  Section  4.  Section  5 presents 
the  results  which  can  be  inferred  from  the  binary  (presence/absenc< ■ 
of  contaminant)  data  from  each  observation  volume  along  the 
sampled  segment.  The  relationship  between  the  theoretical 
considerations  and  the  experimental  results  is  presented  in 
Section  6. 


2.  EXPERIMENTAL  FACILITY  AND  TECHNIQUES 

2.1.  Introduction 

A turbulent  motion  in  the  absence  of  a moan  velocity  was 
generated  by  passing  a bi-planar  grid  through  an  enclosed  volume. 
Light  scattering  techniques  were  used  to  sense  the  presence  or 
absence  of  the  submicron  particles  which  served  as  a marking 
agent  for  the  gas  originally  contained  in  the  lower  half  volume, 
subsequently  referred  to  as  gas  D.  A rapid  scan  of  the  mixing 
region  was  made  possible  by  an  optical-mechanical  system.  Light 
scattering  intensity  measurements  were  sequentially  recorded 
from  short  (0.36mm),  small  diameter  (0.2  mm)  portions  of  a 74  mm 
segment  of  the  incident  laser  beam.  The  physical  apparatus  to 
create  the  mixing  field  and  to  execute  the  measurements  is  de- 
scribed in  detail  in  the  following  sub-sections.  It  should  be 
noted  that  the  rapid  scan  represents  a unique  experimental 
technique.  The  gas  handling  system  that  was  required  for  intro- 
ducing superheated  Freon  12  and  22  into  the  mixing  chamber  is 
described  in  Appendix  A. 

2.2.  Details  of  the  Experimental  Facility 

A schematic  representation  of  the  mixing  chamber  is  shown 
in  Figure  1.  The  mixing  chamber  was  constructed  of  1.25  cm 
plate  stock  of  dimensions  (0. 46x0. 46x0. 61)m;  the  large  I-beams 
provided  a rigid  support  for  the  apparatus.  The  incident  light 
ports  were  covered  with  17.8  cm  diameter  coated  glass  discs; 
the  viewing  port  was  a 12.7  cm  hole  which  accommodated  the 
collection  lens. 

A 0.159-cm  thick  divider  plate  was  located  at  the  midsection 
of  the  chamber  and  separated  the  two  gas  volumes  (A  and  B)  prior 
to  the  initiation  of  the  experiment.  The  plate  extended  through 
the  end  wall  of  the  chamber  which  was  fitted  with  four  Teflon 
guide  surfaces  to  inhibit  the  flexing  of  the  plate  during  ac- 
celeration and  deceleration.  The  grid/splitter  plate  system 
was  supported  on  four  linear  ball  bearings  which  provided  a 
smooth  but  rigid  transport  of  the  plate  and  grid  to  the  forward 
section  of  the  chamber  at  the  initiation  of  the  experiment.  The 
square  mesh,  biplanar,  round  rod  grid  was  constructed  of  aluminum 
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tubing  pinned  at  the  bar  cross-sections;  it  was  permanently 

fixed  at  the  end  of  the  plate.  The  diameter  and  mesh  size  were  ] 

selected  to  provide  a solidity  ratio  of  0.34;  the  tubing  diameter 
was  0.476  cm  with  a mesh  width  of  2.54  cm.  Figure  2 provides 
a photograph  of  the  opened  mixing  chamber.  The  splitter  plate 
has  been  removed;  the  grid  and  the  ends  of  the  four  guide  ways 
are  visible. 

A schematic  representation  of  the  acceleration/deceleration 
system  is  shown  in  Figure  3.  The  quick  acting  solenoid  valve 
allowed  the  rapid  pressurization  of  the  pneumatic  side  of  the 
cylinder.  When  the  unit  is  operated,  the  piston  and  the  attached 
grid  and  splitter  plate  accelerated  until  the  back  pressure  in 
the  hydraulic  fluid  was  sufficiently  large  to  create  a zero  nc 
force  on  the  piston.  The  deceleration  was  accomplished  by  the 
gradual  seating  of  a tapered  plug  attached  to  the  end  of  the 
activating  shaft.  This  facility  provided  a relatively  uniform 
acceleration  from  0 to  6.22  m/sec  in  8.9  cm,  a "constant"  speed 
through  the  chamber,  and  a deceleration  over  the  final  2.5  cm 
of  travel.  This  motion  was  measured  by  attaching  a hot  wire 
(d  = 20pm)  on  the  moving  assembly.  A linearizer  was  used  and 
an  exponent  was  chosen  to  match  the  upper  velocity  range  with 
the  "best"  linear  response.  A plot  of  the  velocity  versus  time 
is  shown  in  Figure  4.  See  Appendix  B for  the  detailed  design 
and  the  computer  simulation  considerations  of  the  acceleration 
deceleration  system. 

2.3.  Optical  System 

A schematic  representation  of  the  optical  system  is  pre- 
sented in  Figure  5.  A Coherent  Radiation  CR-5  Argon  Ion  laser 
was  used  as  the  light  source.  The  laser  was  tuned  for  the 

O 

4880  A line  with  an  output  power  level  of  0.4  watts.  The  2-mm 
diameter  laser  beam  was  focused  with  a 31.75-cm  focal  length 
spherical  lens  and  directed  by  front  surface  mirrors  from  the 
top  to  the  bottom  of  the  chamber  through  two  transparent  windows. 

The  collection  optics  consisted  of  two  aerial  photography 
camera  lenses  (focal  number  of  4.0,  focal  length  of  30.5  cm) 
which  were  placed  in  close  proximity  in  order  to  achieve  a 
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1.36  magnification  of  the  oriqinal  image.  This  represented 
the  maximum  magnification  and  consequently  the  minimum  scan 
length  with  the  available  optics.  These  conditions  were 
necessary  since  the  smallest  scales  of  the  turbulent  motion 
were  estimated  to  be  0.1  mm.  The  resulting  scan  in  the  mixing 
chamber  covered  7.38  cm.  The  first  lens  was  placed  flush  with 
the  inside  surface  of  the  12.7-cm  hole  in  the  side  of  the  chamber. 
This  eliminated  the  need  for  a finite  thickness  window  which 
would  have  caused  spherical  aberrations  of  the  scattered  light. 

The  second  lens  was  coaxial  with  the  first  and  was  moved  axially 
such  that  the  image  of  the  incident  laser  beam  could  be  focused 
onto  the  face  of  a rotating  disk. 

The  high  strength  aluminum  disk  served  as  an  optical-mechanical 

gate  by  only  allowing  the  transmission  of  light  from  a scattering 

volume  defined  by  the  intersection  of  the  recreated  beam  and  a 

slit  cut  in  the  face  of  the  disk.  The  slits  were  machined 

following  the  contour  r = ro  + c (K  where  c^  = 7.14  cm/radian, 

rQ  = 6.99  cm,  0 <_  0 80'  and  0^  = 0 + ^ ( i — 1 ) . The  disk  was 

rotated  at  constant  speed  (4286  rpm)  by  a 1/2  h.p.  motor. 

Since  dO/dt  = constant,  dr/dt  is  also  constant  and  the 

optical-mechanical  gate  provides  a constant  rate  scan  of  the 

recreated  laser  line  over  the  radial  span  r < r < r +10  cm. 

o — — o 

The  scan  is  centered  at  the  midplane  of  the  mixing  field. 

Each  scan  was  completed  in  3.11  msec;  the  scans  were  repeated 
each  3.5  msec.  Figure  6 presents  a detailed  representation 
of  the  disc  and  slits.  Some  difficulty  was  encountered  when 
trying  to  cut  a slit  width  of  0.025  cm  on  the  face  of  the 
disk.  The  machining  operation  left  a groove  which  varied 
linearly  from  0.06  cm  to  0.025  cm.  Following  numerous 
attempts  to  rectify  this  situation,  a solution  was  achieved 
by  placing  Mylar  tape  over  the  machined  groove  and  cutting  a 
slit  with  two  razor  blades  displaced  by  0.025  cm.  One  of  the 
edges  on  the  disk  was  used  as  a guide.  Microphotographs  were 
taken  of  the  slit  variation  after  the  data  had  been  accumulated; 
a sample  is  shown  in  Figure  7.  The  slit  widths  determined 
from  these  photographs  are  presented  in  Table  1.  The  height  of 
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the  scattering  volume  is  a function  of  the  magnification  ratio 

(1.36),  the  diameter  of  the  laser  beam  and  the  intersection 

angle  between  the  slit  and  the  recreated  laser  beam.  Table 

1 also  presents  the  scattering  volume  heiqhts  for  representative 

locations.  The  diameter  of  the  scattering  volume  was  established 

by  the  focusing  lens.  Johnson  11971)  provides  an  analytSyal 

description  of  the  minimum  diameter,  d^,  and  the  length, 

Az^,  over  which  the  diameter  is  maintained.  The  minimum 

diameter  d is  located  at  Az  from  either  end  of  the  "constant" 
o _o 

diameter  region.  2 

A f 

do  = IT  "IT  £ °-1  Inm  (2.3.1) 

O 

where  f = focal  length  of  lens  (31.75  cm),  \,  = 488 OA, 

^ Laser 

d^  = 2 mm. 

Az  = 4(J-)2  A.  = 5 cm  (2.3.2) 

O Cl  _ l_i 

Lj 

The  total  scan  length  was  7.38  cm  and  205  data  points  were  ac- 
quired along  the  scan.  The  distance  between  centers  of  the 
sample  volumes  was  slightly  greater  than  the  nominal  length 
of  the  volumes,  namely  0.36  versus  0.3  mm;  hence  no  known  spat  i 
overlapping  between  samples  was  incurred. 

The  alignment  procedure  was  quite  critical  and  great  care 
was  taken  to  align  the  laser  beam  parallel  to  the  disk  and  the. 
to  focus  a point  source  onto  the  face  of  the  rotating  disk. 

The  latter  was  accomplished  by  adjusting  the  focal  length  of 
the  second  lens  until  the  pulse  width  from  the  photomultipliei 
tube  was  a minimum.  This  alignment  procedure  allowed  less 
than  0.1  degrees  alignment  difference  between  the  incident  and 
focused  laser  beams. 

2.4.  Particle  Generators 

Two  particle  generators  were  developed  in  the  course  of 
this  investigation.  The  first,  which  was  quite  satisfactory 
for  use  with  (at  least)  air  was  not  suitable  for  use  with  Freon. 
Doth  generators  will  be  described  in  this  section. 

The  generator  was  required  to  produce  approximately 
10l>  particles/cc  of  carrier  gas;  this  concentration  is  necessary 
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to  achieve  the  nominal  value  of  20  particles/scattering  volume. 

Further  requirements  include: 

i)  particle  sizes:  0.5*  - 0.3  microns  (imi) 

ii)  non-toxic  particles 

iii)  method  of  generation  should  be  independent  of 
carrier  gas 

iv)  particle  dimensions  should  be  stable  over  the 
lifetime  of  the  experiment. 


The  first  generator  satisfied  requirements  (i)  and  (ii);  it  is 
shown  in  Figure  8. 

The  particles  were  formed  from  soybean  oil  which  was  intro- 
duced as  droplets  into  the  stream  of  carrier  gas.  An  oil  lubri- 
cator, previously  used  in  a compressed  air  system,  was  used  to 
create  the  oil  drops.  The  particle  droplets  for  the  experiment 
were  created  by  the  successive  impacts  of  the  initially  (large) 
drops  on  the  pressed  wire  screens  (of  diminishinq  open  area) 
and  the  subsequent  straining  of  the  liquid  by  the  shear  of  the 
higher  speed  gas.  An  electron  microscope  photograph  of  the 
resulting  droplets  (see  Figure  9)  indicated  that  droplet 
diameters  of  0.2  - 0.4  microns  were  achieved  with  this  generator. 
Unfortunately,  the  droplets  so  produced  in  Freon  appeared  to 
evaporate  (or  at  least  to  change  their  surface  charge  character- 
istics) . No  particle  scattering  was  observed  when  this  scheme 
was  used  with  Freon  as  the  carrier  gas. 

The  second  particle  generator  satisfied  all  four  criteria; 
it  is  shown  schematically  in  Figure  10.  This  generator  makes 
use  of  the  Rayleigh  instability  of  a capillary  jet  and  the 
subsequent  formation  of  droplets  to  create  the  desired  marking 
particles.  The  stability  of  the  particle  was  ensured  by  using 
a liquid  solution  of  10  percent  acetic  acid*  in  distilled  water. 
When  the  liquid  water  evaporated  following  the  droplet  formation 
a crystalline  residue  was  left  to  serve  as  the  scattering  par- 
ticle. It  was  inferred  that  the  0.45-um  micro  porous  filter 
provided  a capillary  jet  of  the  same  diameter  at  its  downstream 
face.  A capillary  jet  is  unstable;  the  analysis  of  Lord  Rayleigh 
(1892]  was  used  to  infer  that  0.9-um  droplets  were  formed.  An 

*Acetic  acid:  CH,  COOH,  m.w.  60.03 
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unambiguous  measure  of  the  particle  sizes  was  not  obtained.  An 

★ 

electron  microscope  photograph  of  the  particles  which  settled 
from  the  mixing  chamber  is  shown  in  Figure  11.  A long  time 
(24  hr.)  was  required  for  the  settlinq  and  some  agglomeration 
of  the  initially  distinct  particles  is  suspected.  This  tech- 
nique is  vulnerable  to  oversize  droplets  which  impinge  upon 
the  upstream  face  of  the  filter;  the  apparent  difficulty  is  the 
merging  of  the  capillary  jets  and  hence  the  formation  of  droplets 
of  an  excessive  size.  The  generator  system  appeared  to  operate 
the  same  for  the  different  carrier  gasses ....  Freon  12  and  22 
and  Air.  An  improvement  in  the  efficiency  of  the  system  would 
be  desirable  since  an  asymptotic  state  is  reached  when  the 
particle  contaminated  gas  is  recycled  through  the  compressor 
and  particle  generator. 

2.5.  Signal  Processing 

The  scattered  light,  which  is  collected  and  focused  onto 
the  photomultiplier  tube  by  the  optical  system,  results  in  a 
voltage  signal  which  can  be  used  to  signify  the  presence  or 
absence  of  the  marking  particles.  The  characteristics  of  the 
photomultiplier  tube  and  the  signal  processing  electronics  are 
described  in  this  section. 

An  RCA  4518  photomultiplier  tube  was  used  for  the  experiments. 
2 

This  tube  has  a 20  cm  photosensitive  area  and  a quantum  efficiency 
of  18  percent.  The  10K  ohm  load  resistor  and  the  associated 
capacitance  of  the  output  circuit  resulted  in  a 5-usec  time  con- 
stant. The  tube  was  operated  with  900  volts  across  the  dynode 
chain;  the  corresponding  cathode  sensitivity  and  current  ampli- 
fication were  0.06  amps/watt  and  8 x 10J  respectively. 

The  particulate  seeding  was  nominally  controlled  to  pro- 
vide 20  particles  per  scattering  volume  for  the  marked  fluid. 

With  this  seeding  density  and  for  an  incident  laser  power  of 
0.4  watts,  the  optical/electronics  system  delivered  a current 
of  approximately  one  micro  amp.  This  current  level  corresponds 

4 

to  approximately  10  photoelectrons  for  each  5-ysec  sample  period. 
This  relatively  large  number  guarantees  that  the  photomultiplier 

* A 1%  solution  of  Fhototungstic  acid  (PTA)  was  added  to  the 
distilled  water  to  facilitate  the  preparation  of  the  slide. 


tube  and  associated  circuitry  will  not  be  adversely  influenced 
by  the  dark  current  and  that  a statistically  reliable  represen- 
tation of  the  particles’  presence  can  be  expected. 

The  photomultiplier  tube  signal  was  amplified  IPAR  101)  by 

2 

10  and  input  to  an  electronic  comparator  circuit.  The  com- 
parator level  of  0.06  volts  was  sufficiently  above  the  "system 
noise/Rayleigh  (molecular)  Scattering  dark  current"  level  of 
0.02  volts  such  that  no  spurious  readings  were  attributed  to  these 
effects.  Voltage  levels  of  greater  than  O.On  volts  were  con- 
sidered to  represent  the  presence  of  particulate  matter  and 
hence  the  presence  of  gas  B in  the  scattering  volume.  A minimum 
value  of  approximately  0.2  volts  was  observed  for  the  region  of 
particulate  matter  when  the  experiment  was  initiated;  hence, 
the  signal-to-noise  level  was  approximately  10  and  the  minimum 
signal  level  exceeded  the  threshold  by  approximately  3.  The 
voltage  amplitude  associated  with  the  presence  of  the  particulate 
matter  did  not  exhibit  a unique  value;  the  amplitude  fluctua- 
tions were  inferred  to  be  approximately  *50  percent  of  the 
average  amplitude  for  those  scattering  volumes  which  contained 
the  marked  fluid.  Consequently,  it  was  not  possible  to  set  a 
criterion  for  the  discrimination  of  a partially  filled  scattering 
volume.  There  are  two  principal  factors  which  contribute  to 
the  amplitude  fluctuation  when  the  scattering  volume  is  com- 
pletely filled  with  marked  fluid:  (i)  the  nominal  seeding  of 
20  particles  per  scattering  volume  results  in  a standard  devia- 
tion of  *4.5  particles  or  *23  percent  (assuming  a Poisson 
distribution  for  the  number  of  particles  in  a scattering  volume) 
and  (ii)  the  variation  in  the  scittering  volume  dimensions  over 
the  four  scans  was  approximately  *25  percent. 

Figure  12  presents  an  oscilloscope  trace  of  the  photo- 
multiplier tube  output  immediately  following  the  passage  of 
the  grid  and  splitter  plate.  The  relevant  voltage  levels  are 
noted  on  this  trace.  Figure  13  shows  a block  diagram  of  the 
data  processing  circuitry  which  achieves  the  analog  to  digital 
conversion  and  which  synchronizes  the  data  from  the  P.M.  tube 
and  rotating  disk.  The  comparator  curcuit  which  was  used  to 
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process  the  analog  signal  from  the  photomultiplier  tube  is 
shown  in  Figure  14. 

Each  scan  of  the  mixing  field  is  initiated  by  the  strobe 
signal  from  the  phototransistor.  The  phototransistor  output 
is  activated  by  the  incident  HeNe  (low  power)  laser  beam 
passage  through  the  slit.  This  scheme  for  initiating  the  data 
scan  is  estimated  to  provide  repeatability  of  0.05  mm  (or  14 
percent  of  one  scattering  volume)  in  the  spatial  locations  for 
the  separate  traverses.  Figure  15  shows  the  binary  output  of 
the  mixing  field  for  one  experimental  run. 


3.  ANALYTICAL  AND  BACKGROUND  CONSIDERATIONS 
3.1  Description  of  the  Problem 

It  is  useful  to  consider  the  present  study  in  terms  of  an 
idealized  experiment  and  the  differences  between  it  and  the 
actual  experiments  which  were  executed.  The  idealized  experiment 
involves  a homogeneous,  isotropic  turbulence  field  (of  indefinite 
extent)  at  the  initiation  of  the  experiment  (t=0) . The  initial 
concentration  field  is  specified  by  the  conditions 


r(x,y,z,0)  = 1 for  z < 0 
r(x,y,z,0)  = 0 for  z > 0 


(3.1.1) 


where  T = 1 indicates  the  presence  of  the  non-dif fusive  contaminant 
in  the  differential  sized  volume  located  at  the  (x,y,z)  point. 

The  subsequent  r(x,y,z,t)  distributions  result  from  the  convective 
transport  of  the  turbulence  field.  The  actual  experiment  is  rather 
different  in  certain  of  its  detailed  features.  The  passage  of 
the  splitter  plate  trailing  edge  provides  an  extremely  rapid 
mixing  of  the  upper  and  lower  fluid  regions  over  a scale  several 
times  the  plate  thickness  (t  ^afce  = 1.56  mm).  This  provides  an 
actual  "initial  condition"  which  is  a well-mixed  region  for  ±7mm. 
Detailed  considerations  of  the  initial  condition  are  presented 
in  Section  5.  The  balance  of  the  initial  turbulence  field  is 
that  which  results  from  the  passage  of  the  biplanar  grid. 

The  x-component,  ensemble  mean  velocity  will  be  positive 
(i.e.,  the  fluid  motion  will  be  toward  the  receding  grid)  in 
the  wake  regions  of  the  grid  and  splitter  plate.  The  velocity  u 
in  the  open  regions  will  be  negative.  The  closed  mixing  chamber 
results  in  the  integral  constraint: 

/ / u dy  dz  = 0 (3.1.2) 

for  all  y and  z in  the  interior  of  the  chamber. 

The  differential  forms  of  the  mass  and  momentum  conservation 
equations  are  instructive  in  the  interpretation  of  the  initial 
region  behavior.  Since  the  subject  gases  behave  as  incompressible, 
Newtonian  fluids,  these  equations  are  written  in  their  ensemble 
averaged  form  as 
9u , 

(3.1.3) 
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(3.1.4) 


The  turbulent  stresses  of  ( 3 . 1_.  4)  will  serve  to  accelerate 
(■^  > 0)  and  to  decelerate  (-|^  < 0)  the  fluid  behind  the  grid 
members  and  in  the  open  regions  respectively.  The  natural  symmetry 
of  the  grid  and  plate  leads  to  planes  of  zero  v and  w and  the 
conservation  of  mass  equation  can  be  used  to  infer  non-zero 
v and  w values  between  these  planes  in  response  to  the  observed 
3u/3x  values.  For  example,  for  the  y = 0 plane  between  two 
grid  members  (i.e.,  v = 0 for  all  z) , w is  zero  at  the  mid-plane  of 
a horizontal  rod,  it  is  negative  in  the  region  between  the  plane 
of  the  rod  and  the  center  plane,  and  it  is  zero  at  the  midplane 
between  the  rods.  The  presence  of  this  mean  motion  following 
the  passage  of  the  grid  suggests  that  convective  transport 
effects  are  present  in  the  actual  experiment  that  would  not  be 
present  if  the  idealized  initial  condition  of  homogeneous  tur- 
bulence were  obtained.  The  physical  extent  of  the  inhomogeneous 
region  can  be  estimated  from  the  experimental  results  of  Sato 
and  Yamamoto  [1974].  Their  hot-wire  traverses  downstream  of 
a biplanar,  round  rod  grid  of  solidity  0.44  (vs.  0.34  for  the 
present  study)  indicate  that  homogeneity  is  obtained  in  approxi- 
mately 3 mesh  lengths.  This  corresponds  to  approximately  the 
fourth  scan  for  the  present  study. 


The  turbulence  kinetic  energy  equation  with  q = 
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(3.1.5) 


makes  explicit  that  the  turbulence  kinetic  energy  is  generated 
in  this  initial  inhomogeneous  region  and  decays  (balance  of  first 
and  last  terms)  as  a result  of  dissipation  effects  thereafter. 

The  maximum  magnitude  of  the  turbulence  kinetic  energy  can  be 
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established  in  terms  of  the  work  done  by  the  grid/splitter 

plate  movement  through  the  chamber.  Let  C'D  be  the  drag  coefficient 

for  the  physical  structure  used  to  create  the  turbulence  field. 

The  total  work  done  is 
W =*  / ? • 3s 

C V 2A  -♦  -+ 

Since  F = Dp  o G and  ds  = V dt 
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here  t is  the  time  required  to  move  the  assembly  across  the 
chamber.  Since  the  grid  movement  is  controlled  by  the  hydraulic 
actuator  and  since  this  motion  is  independent  of  the  drag  force 
acting  on  the  grid,  the  initial  kinetic  energy  states  for  the 
air  and  Freon  gases  will  be  in  the  ratio  of  their  CD  values. 
Assuming  that  the  CD  value  is  set  by  the  grid  elements,  the 
data  of  Friehe  and  Schwarz  [1970)  can  be  used  to  establish  the 
ratio 
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It  is  well  known  that  in  a parallel  flow  5 is  the  turbulence 
velocity  component  which  is  first  and  predominantly  energized  by 
the  destruction  of  the  mean  flow  energy;  see  Corrsin  [1957).  It 
is  expected  that  this  anisotropic  condition  is  present  in  the 
current  experiment.  A very  considerable  data  base  is  available 
for  the  flow  behind  biplanar  grids.  Comte-Bellot  and  Corrsin 
[1966]  have  summarized  these  results  including  those  of  Wyatt  [1955) 
which  are  for  a configuration  very  similar  to  that  of  the  present 
study. 

Following  the  establishment  of  a homogeneous  turbulence 
field  and  prior  to  the  final  period  of  decay,  the  available 
experimental  data  suggest  that  a power  law  solution  to  (3.1.5)  is 
appropriate.  Comte-Bellot  and  Corrsin  [1966]  organize  the  ex- 
perimental results  in  terms  of 


U 


o 


2 


. U t U t 
A r O O Q| 

1 M M ' 


n 


M 


(3.1.7) 
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which  includes  the  presence  of  an  apparent  origin  to  account 
for  the  initial  period  of  inhomogeneous  turbulence.  Using  the 
data  of  Wyatt  [1955]  a graphical  representation  of  the  pertinent 
turbulence  quantities  in  the  present  experiment  has  been  pre- 
pared, see  Figure  16. 

Data  for  round  rod,  biplanar,  grids  show  that  the  apparent 
origin  (tQ)  shifts  to  smaller  values  as  the  velocity  is  increased 
and  mesh  size  is  kept  constant;  see  Comte-Bellot  and  Corrsin 
[1966].  Wyatt's  data  suggests  that  Uoto  decreases  from  5 to  4 

4 M 4 

for  a corresponding  range  of  10*  to  6x10  . These  are  the 

relevant  RM  for  the  present  experiment.  However,  it  is  clear, 

from  other  data  reported  by  Comte-Bellot  and  Corrsin,  that 

well  defined  values  of  Uofco  do  not  exist. 

M 

In  their  study  of  the  decay  relations  behind  a biplanar 

grid,  Comte-Bellot  and  Corrsin  have  reported  that  there  was 

H H 

no  apparent  change  in  the  decay  laws  for  ^ = 9 and  ^ = 18, 
where  H represents  the  lateral  dimension  of  the  wind  tunnel. 

The  corresponding  value  for  this  experiment  is  — = 18.  Uzkan 
and  Reynolds  [1967]  have  reported  an  experiment  in  which  the 
decaying  turbulence  field  was  bounded  (on  one  side)  by  a wall 
which  moved  at  the  free  stream  velocity  UQ.  The  contamination 
of  the  interior  turbulence  by  the  wall  effects  is  seen  to  pro- 
ceed at  a reduced  rate  for  this  condition.  Their  data  show 
that  it  would  require  54  seconds  for  the  boundary  influence 

4 

to  migrate  5 cm  into  the  free  stream  at  Rj^  = 10  . Since  times 
of  less  than  0.8  seconds  are  of  interest  in  the  present  ex- 
periment, it  is  inferred  that  there  is  no  influence  of  the 
boundaries  on  the  measurements  reported  herein.  Bradshaw  (1977) 
has  reported  that  experiments  of  this  type  are  sensitive  to  the 
Reynolds  number;  however,  the  influence  of  the  stationary  wall 
will  not  influence  either  the  Freon  or  the  air  results. 


3.2.  Scalar  Concentration  Field 

3.2.1.  The  Ensemble  Average  Concentration  Distribution 

The  transport  equation  for  a scaler  contaminant,  C,  is 


C 


a2c 

3x  j 3x  j 


(3.2.1) 


where  C is  presumed  to  follow  a Fickian  law  of  diffusion.  For 
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the  present  discussion,  it  is  useful  to  consider  the  idealized 
condition  in  which  D - 0;  this  is  closely  approximated  in  the 
experiment  as  shown  by  the  analyses  in  Section  4. 

For  the  condition  that  D = 0,  the  mean  concentration  of 
the  nondiffusive  contaminant  F(i.e.,  is  given  by 
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(3.2.2) 


Since  <u^>  0,  except  for  the  initial  period  of  inhomogeneity, 

the  expression  for  this  restricted  condition  is  of  interest. 
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(3.2.3) 


In  3.2.3  the  turbulent  transport  term  is  written  as  a 
product  of  the  turbulent  diffusivity  and  the  gradient  of  the 
mean  concentration  field.  This  is  based  on  an  ad  hoc  analogy 
to  molecular  transport  modeling,  Corrsin  [1974] . 

For  the  problem  under  consideration,  namely,  homogeneous 
turbulence,  D ( can  be  assumed  to  be  independent  of  the  direc- 
tion of  diffusion;  consequently,  (3.2.3)  can  be  written  as 


a<r> 

at 


n,  (t) 


with  boundary  cond i t i ons 


( 3 . 2 . 4 . a ) 


vf  ( z , 0 ) s = 0 , z:-0 
-f  (z,0)>  = 1 , Z'O 


(3.  2.4.b) 


The  general  solution  of  this  equation  in  a semi-infinite 
region  Zv«-)  is  given  by  Ozisik  119681  . 

v 1 ( z , t ) i ^ I exp  -{Z~*  ^ 1’  ( z 1 0 ) d z ' 

{ 4 n / D ( t * )dt ' } {4  / D (t)dt  ' } 

o o 1 

(3.2.5) 

where  I(z'O)  is  the  initial  concentration  distribution. 

^ Note  that  both  ( ) and  «•  ( ) s will  be  used  to  designate  an 
ensemble  average  of  ( ) . 
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The  same  result  for  <T(z,t)>  can  be  obtained  by  an  inde- 
pendent analysis  which  is  due  to  Batchelor  [1956] . This 
alternative  development  is  considered  in  order  to  clarify  the 
assumptions  and  restrictions  which  are  used  in  each  analysis. 

For  the  conditions  of  homogeneous  turbulence,  without 
mean  velocity  gradients,  it  is  reasonable  to  assume  that  the 
three-dimensional  particle  displacements  follow  a Gaussian 

distribution  with  standard  deviations  (x  , y , z ) . The 

P P P 

probability  density  function  of  the  particle  displacements 
is  therefore 


Q(x,y , z, t) 


( 2 it  ) 5$ 


7 ^2 2~Ij  exP 
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For  a scalar  contaminant  which  indelibly  marks  the  fluid 
element  and  for  one  whose  presence  is  identified  by  r = 1 
(and  absence  by  T = 0) , the  ensemble  average  concentration 
distribution  of  the  contaminant  is  directly  related  to  the 
displacement  of  the  original  particle  distribution  as 


CO 

<T(x,y,z,t)>  = ///  Q (x-x ' ,y-y ' , z-z ' , t)  T (x ' ,y ' , z ' , 0 ) dx'dy'dz' 

— CO 


(3.2.7) 

where  T (x ' , y ' , z ' , 0 ) is  the  initial  concentration  distribution. 

For  the  experiment  under  consideration,  the  initial  distribu- 
tion of  the  contaminant  occupies  the  volume  (x,  y,  zN0) ; 
therefore,  <T(x,y,z,t)>  = <T(z,t)> 
and 


<T  ( z , t)  > = / 1 v exp 

-«>  ( 2ttz 
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(z-z'  ) 
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(z*  0) 


dz  ' 


(3.2.8) 


This  result  suggests  that  the  mean  square  displacement 
of  a fluid  particle,  which  represents  the  dispersive  property 
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of  the  turbulence  field,  governs  the  behavior  of  the  ensemble 
average  concentration  distribution.  The  mean  square  displace- 
ment can  be  developed  using  Lagrangian  considerations;  this 
development  is  presented  in  Appendix  D. 

The  above  equation  is  similar  to  (3.2.5),  where  by  iden- 
tity. 
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2 / 
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Dt  (t' )dt* 


(3. 2. 9. a) 


or 


DT(t) 


( 3 . 2 . 9 . b) 


This  result  was  observed  by  Hinze  [1975].  Hence  the  turbulent 
diffusivity  can  be  calculated  from  the  time  rate  of  change  of 
the  particle  mean  square  displacement. 

The  mean  concentration  equation,  (3.2.8),  can  be  written 
in  terms  of  the  error  function  by  defining  new  independent 
variables;  see  Ozisik  [1968] . Let 

« = — Z~o  i,  and  dz'  = (2  z 2)^  da  (3.2.10) 

( 2z  P 

P 

Utilizing  the  boundary  conditions  from  (3. 2. 4. a),  the  follow- 
ing expression  is  obtained  for  the  ensemble  average  distribu- 
tion; viz.  , 

<T(z,t)>  = \ (1  + erf  (z) ) (3. 2. 11)  + 

For  the  experimental  data  under  consideration  a Gaussian 
cumulative  distribution  function  which  is  identical  to  (3.2.8) 
will  be  used  to  describe  T(z,t),  namely, 

z * 2 

<r(z,t)>  = ^ — / exp  - j.a~<Z>)  da  (3.2.12) 

( 2 tt  ) o 2 a 2 
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The  separate  symbols;  o and  z , will  be  maintained  to 

m p 


+The  error  function  of  (3.2.11)  is  defined  as 
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respectively  represent  the  best  fit  of  a Gaussian  C.D.F.+  to 
the  experimental  data  for  <T(z,t)>  and  the  standard  deviation 
of  the  particle  displacements  in  (the  assumed  Gaussian  dis- 
tribution of)  the  flow  field. 

The  rate  of  transport  of  the  marking  contaminant  across 
the  midplane  of  the  mixing  chamber  can  be  expressed  in  terms 
of  the  variables  used  in  the  above  considerations.  The  ensemble 
average  flux  term  of  the  contaminant  across  the  midplane  is 
<wy>z=0  and,  by  using  (3.2.3),  (3.2.5)  with  the  r(z',o)  given 
by  (3.2.4.b),  and  (3.2.9.b),  the  expression  for  <w>>z=o  becomes 


<wy  > 


z=o 
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(2tt)  ^ 


d.ZE 

dt 


(3.2.13) 


The  total  flux  of  <T>  during  the  elapsed  time  period  o£t'£t 
may  be  expressed  as  the  integral  of  (3.2.13);  viz.. 
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> dt' 
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(3.2.14) 


3.2.2.  The  Ensemble  Average  Fluctuation  Intensity  Distribution 


is 


The  mean  square  concentration  fluctuation,  <c  >,  equation 
3 2 9<c>  92  „ 9c  9c  „ 32c2 

st  <c  >=  - <ujc  Tnir~  " S3r  <c  uj>  * D <j^r  7x7>  " D 

j 3 3 3 3 3 


(3.2.15) 


where  c = C - C in  which  the  turbulent  effects  serve  to  pro- 
duce and  convectively  transport  the  concentration  fluctuations 

(i.e.,  the  first  two  terms  on  the  right  side  of  (3.2.15))  and 

2 

the  molecular  effects  destroy  and  diffuse  the  <c  > by  the 
action  of  the  second  two  terms.  In  the  absence  of  diffusion, 
i.e.,  for  D = 0,  a balance  is  possible  if  the  production  and 
transport  are  in  balance.  This  condition  occurs  at  the  mid- 
plane of  the  mixing  region  as  shown  by  the  following  analysis. 

If  D = 0,  then  the  probability  density  function  for  r, 
Pp(n,  t) , can  be  represented  by  two  Dirac  functions  at  a = 0 


t Cumulative  Distribution  Function 
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and  a 


1 


such  that  the  mean  value  of  r is  obtained  as 
1 

<r>=/  Pp  (a,  t)ada»Pr(lf  t)  ( 3 . 2 . 1 e» ) 

o 


The  mean  square  fluctuation  is 

1 o 

<T>  = / (a-<T>)^  Pr  (a , t ) d a 

o 

- <r>2  Pr(o.t)  - n -<r>) 2 Pra.t)  {3.2.17) 

= <r>(i-<r>) 

Hence,  at  the  centerline,  <y  > is  constant  for  all  time  (since 
<r>  = 0.5);  its  value  is  0.25. 


2 

Using  the  results  of  (3.2.17),  <y  > can  be  written  as 

<Y2>  = 0.25  {1  - erf2  (3.2.18) 

( 2z  ) 

P 

Measurements  by  Foss,  Schlien  and  Corrsin  [1972]  in  the  mixing 
field  downstream  of  a half-heated  grid  showed  that  the  above 
estimate  for  the  temperature  fluctuation  was  250%  larger  than 
the  experimental  data.  Hence,  the  molecular  diffusivity  plays 
an  important  role  in  such  flows  even  though  the  Gaussian  C.D.F. 
of  equation  (3.2.8)  was  closely  approximated  by  the  experimental 
data  for  the  mean  temperature  distribution.  Analytical 
descriptions  of  the  fluctuation  intensity  for  the  diffusive 
condition  have  been  developed  by  Libby  [1975]  and  the  above 
investigators . 

3.2.3.  Dispersion  Considerations 

The  T (z,t)  distributions  which  have  been  recorded  in  the 
present  work  represent  the  effects  of  turbulent  transport.  It 
is  interesting  to  note  that  the  non-dif f usive  nature  of  the 
sensed  contaminant  results  in  a direct  correspondence  between 
<T(z,t)>  and  the  turbulent  transport.  This  is  in  contrast 
with  studies  which  use  thermal  tracers  as  the  marking  element. 
Taylor  [1935]  was  the  first  to  analytically  describe  such 
transport  by  the  application  of  Lagrangian  concepts.  Using 
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as  the  velocity  of  a fluid  dynamic  particle,  the  Lagrangian 
auto  correlation  for  stationary  turbulence  can  be  expressed  as 

<v  • (t) y. (t+a) > 

— j (3.2.19) 

<y  > 

and  the  mean  square  displacement  of  a fluid  element  from  its 
point  of  origin  in  a stationary,  homogeneous  turbulence  field  is 

? o t* 

<xi  > = 2<yi  > / dt'  / R (a)da  (3.2.20) 

o o 

for  stationary  turbulence,  the  limiting  cases  of  short  time: 

2 2 2 

t+0,RL(a)  = 1 and  <x  > = 2<y  > t and  of  long  time:  t-*-°°, 

CO 

2 2 

/ R (a)da  = A and  <x  > = 2 <y  >At,  are  of  interest. 

A turbulence  field  without  mean  shear  will,  of  course, 
decay.  If  a tagged  particle  is  released  in  such  a field,  its 
dispersion  as  a function  of  time  will  be  influenced  by  the 
decay  of  the  turbulence  kinetic  energy.  A theoretical  descrip- 
tion of  this  physically  realistic  condition  would  require  that 
the  Lagrangian  correlation  function  would  specifically  account 
for  the  decaying  turbulence  kinetic  energy. 

Falco  [1972]  made  direct  optical  measurements  of  the 
velocity  of  small  particles  in  the  decaying  turbulence  field 
behind  a biplanar  grid  in  an  enclosed  water  tank  facility. 

From  these  measurements  he  obtained  a (spherical)  Lagrangian 

correlation  at  Uofc  = 89  and  156.  It  is  apparent  from  these 
M 

measurements  that  there  is  a temporal  growth  of  the  Lagrangian 

scales  as  Uofc  increases;  however,  due  to  the  sparse  nature 
M 

of  the  data  (8  points  for  the  correlation)  it  would  be  dif- 
ficult to  determine  the  quantitative  features  of  this  growth. 
Several  authors,  including  Batchelor  and  Townsend  [1956]  and 
Shlien  and  Corrsin  [1974],  have  formulated  analyses  to  account 
for  the  decay  of  the  turbulence  energy  and  its  relationship 
to  the  dispersion  of  a contaminant.  A new  version  of  such 


The  common  feature 


an  analysis  is  presented  in  Appendix  D.  The  common  feature 

of  these  analyses  is  the  assumption  of  a self  preserving  state 

of  the  Lagrangian  correlation  function,  R (i),  by  the  use  of 
. t L 

a stretched  time  coordinate,  t.  The  equivalent  results,  for 

2 

the  mean  square  displacement,  , of  a fluid  particle  for 

the  short  and  long  time  special  cases,  are  given  below.  The 

dispersion  measure  in  the  z direction  of  the  present  experi- 
2 

ment,  z , has  been  made  nondimensional  using  the  scale  length 

L where 
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L 2 = <v2  (t  - t ) > (t  - t ) 
p l o r o 


(3.2.21) 
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tf  is  the  "release  time"  for  the  tagged  particle,  tQ  is  the 
apparent  temporal  origin  of  the  decaying  turbulence  and  ” 
is  the  Lagrangian  velocity  of  the  tagged  particle.  The  short 
time  result  is,  for 
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(3.2.22) 


Note,  tj  and  t*  are  defined  below;  (3.2.22)  is  also  equation 
D.2.3.  The  long  time  result  is 


+ This  assumption  was  introduced  by  Batchelor  and  Townsend 
[1956] . 
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(3.2.23) 


The  symbol  n represents  the  decay  exponent  of  the  turbulence 
kinetic  energy;  see  equation  (3.1.7).  The  variables  t*  and 

TL  A 

A*  are  defined  as  t*  = pr—  and  A*  = where  tt  and  A are  the 

1 C2  L 

Lagrangian  micro  and  macroscales  in  the  stationary  reference 
frame.  Theoretical  curves  of  3.2.22  and  3.2.23  for  various 
i and  A values  are  presented  and  discussed  in  Article  6. 

Li 

The  constants  and  are  related  to  the  temporal  trans- 
formation in  which  the  Lagrangian  velocity  correlation  (3.2.19) 

g 

is  stationary;  i.e.f  -y^-(RL  (t , t ) ) = 0.  The  temporal  trans- 
formation which  yields  a stationary  condition  for  Rl(t)  for 
r £t  j is 


— = exp  [i — ] 
-r 


(3.2.24) 


where  = (l-^)C^.  The  similarly  defined  temporal  transfor- 


mation for  t s t . is 
3 


= exp  [pr-1 


(3.2.25) 


where  C j.  = (1--})CV  The  symbols  At,  Atr  and  At^  are  the 

times  referenced  to  the  temporal  origin  as  At  = t-tQ, 

At  = t -t  , and  At.  = t.-t  . The  symbol  t.  represents  t 
rro  DDO  1 


The  symbol  represents  the 


25 


juncture  point  between  the  short  and  long  time  formulations. 

Figure  16b  shows  a graphical  representation  of  these  concepts. 

The  following  relationship  defines  the  juncture  point  in  the 

reference  frame: 

At 

r i-n 

„ 2 , 2. 

1-  j - At . At.  ” A * 

1 - (-£)  In  (^)  ~ (^)  = 0 (3.2.26) 

r r 

The  Lagrangian  macroscale,  A,  is  a property  of  the  turbulence 

field  and  is  clearly  related  to  the  dispersive  character  of 

a turbulent  flow.  Previous  experimental  efforts  to  determine 

the  dispersion  have  made  use  of  a passive  scalar  contaminant 

such  as  a low  overheat  thermal  marker.  In  these  studies, 

temperature  tagging  wires  mark  a spot  or  line  of  the  fluid 

and  downstream  temperature  measurements  define  the  mean 

temperature  as  a function  of  the  transverse  distance.  The 
2 . 

value  Zp  is  inferred  from  the  latter.  Such  studies  have  been 
performed  by  Uberoi  and  Corrsin  11953) , Townsend  11956] , and 
Shlien  and  Corrsin  [1974].  The  last  authors  provide  a well 
developed  scheme  for  interpreting  the  influence  of  the  de- 
caying non-stationary  condition  of  the  turbulence.  The  in- 
fluence of  the  non-zero  diffusivity  of  the  thermal  energy  in 
such  experiments  has  been  theoretically  examined  by  Saffman 
[1961].  Saffman's  analysis,  which  is  limited  to  small  times 
from  the  initiation  of  the  mixing,  provides  a corrective 
interpretation  to  the  earlier  Townsend  [1956]  analysis  which 
had  served  a similar  purpose  for  the  original  work  of  Taylor 
[1935].  Chevray  [1977]  has  developed  an  alternate  analytical 
structure  which  extends  the  prior  results  of  Saffman.  Quali- 
tatively, the  present  investigation  is  as  if  the  marking  con- 
taminant were  impressed  upon  a region,  instead  of  a line  or 
spot,  of  the  host  fluid. 

3.2.4.  Fourier  Transform  Considerations 

A useful  description  of  the  scalar  concentration  field 
can  be  gained  by  the  use  of  Fourier  transforms.  Let  the 
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function  B(k)  represent  the  amplitude  of  the  wave  which  is 
related  to  the  spatial  distribution  of  T(x)  as 
-*  00  >►  -*• 

B(k)  = (211)  ///  r (x)  exp  (-i  k*x)  dx  (3.2.27) 

— oo 

Similarly,  the  velocity  component  un  can  be  transformed  to  a 
wave  description  as 

— ► 

-*  oo  -►  -►  -> 

Fn(k)  = (211)  ///  un(x)  exp(-ik-x)  d x.  (3.2.28) 

— OO 

Operating  on  the  conservation  equation  for  T results  in 
(see,  eg.  Yeh  and  Van  Atta  [1973]) 

f^[i4>(k)]  = -i  kj  ///  Tn  (k,  k ’ ) dk’  - D k^0n  (k)  (3.2.29) 

3 

where  4>  = < BB*  > Ak  which  represents  the  intensity  of  the 
concentration  fluctuation  at  wave  numbers  k^  ± Ak^/2.  (B*  is 

the  complex  conjugate  of  B. ) The  function  Tn(k,  k ' ) represents 
the  transfer  of  the  fluctuation  intensity  to  higher  wave  numbers; 
formally,  it  can  be  expressed  as: 

Tn(k)  = U(Ak)3  ki  <F*(k  - k')B*(k')  B(k)  - Fn  (k-k')  B*  (k)  B(k')  > 

(3.2.30) 

This  transfer  term  is  responsible  for  enriching  or  transfer- 
ring energy  (<j>(k))  to  large  values  of  k.  Its  effect  is  limited 
by  the  presence  of  the  last  term  in  (3.2.11)  in  which  the 
concentration  fluctuations  are  "smoothed"  by  the  diffusive 
action  of  the  molecular  diffusivity.  For  the  idealized  ex- 
periment involving  D = 0,  the  striations  of  the  mixed  regions 
will  become  of  finer  and  finer  scale.  Batchelor  [1959]  has 
developed  a theoretical  description  of  the  scalar  field  for 
the  condition  of  small  conductivity  (diffusivity) . His  work 
shows  that  the  <j>  (k)  spectra  falls  off  at  the  lesser  slope  of 
k ^ beyond  the  break  where  Fn(k)  changes  from  a -5/3  to  a 
region  of  more  rapid  decay.  Gibson  and  Schwarz  [1963] , Nye 
and  Brodkey  [1967],  and  Grant  & Nisbet  [1957]  have  observed  the  k_1 
behavior  in  : i)  salt  concentration  and  temperature  measure- 
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ments,  ii)  dye,  and  iii)  temperature  measurements  respectively. 
Interestingly,  there  was  no  significant  difference  between 
the  salt  and  heated  water  results  in  the  Gibson  and  Schwarz 
experiment. 
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4.  DYNAMIC  BEHAVIOR  OF  THE  MARKING  PARTICLES 

The  previous  sections  have  identified  the  basic  premise 
of  the  experimental  technique:  that  a marking  contaminant  may 
be  found  which  faithfully  moves  with  the  same  continuum  fluid 
element  from  the  initiation  to  the  termination  of  the  experiment. 
The  contaminant  particle  selected  must  be  sufficiently  large 
that  it  is  not  displaced  from  the  continuum  element  by  the 
random  agitation  of  the  molecular  motion  and  it  must  be  suffi- 
ciently small  that  it  does  not:  i)  follow  an  inertia  dominated 
trajectory  and  thereby  define  a pathline  which  is  different 
from  that  of  the  continuum  element  of  its  origin,  ii)  influence 
the  density  or  stability  properties  of  the  fluid,  or  iii)  fall 
away  from  the  continuum  element  under  the  action  of  the  gravi- 
tational force. 

The  following  subsections  consider  these  effects  based 
upon  air  as  the  host  fluid  for  submicron  particles  (a  = 0.25pm) 
of  specific  gravity  1.0. 


4.1.  Aerodynamic  Drag 

The  equation  which  describes  the  motion  of  the  fluid 
element  and  the  corresponding  response  of  the  particle  has  been 
analyzed  by  Tchen  [1947] , Hinze  [1975] , Corrsin  and  Lumley  [1956] 
and  Soo  [1967] . The  following  equation  was  extracted  from  Khosla 
and  Lederman  [1973] . 
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(4.1.1) 


where  u^  is  the  vector  velocity  of  the  particle,  u^  is  the 

velocity  of  the  gas,  v the  coefficient  of  kinematic  viscosity, 

M the  coefficient  of  viscosity,  p the  density  of  the  solid 

hr 

*Becker  [1977]  identified  this  technique  as  "marker  nephelometry " 
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material  and  a is  the  radius  of  the  particle.  Tchen  assumes 
that  the  particle  radius  is  small  and 


^ >>  1 (4.1.2) 

2 3 u . 
a l 

3x  . 3x  . 

D 3 

The  most  restrictive  assumption  made  by  Tchen  is  that  during 
the  motion  of  the  particle,  the  same  fluid  element  remains  in 
its  neighborhood.  Thus  he  postulates 


d d 
dtp  - dt  * 


(4.1.3) 


This  assumption  is  required  for  a solution  to  be  realized  so 
that  one  can  achieve  an  estimate  of  the  phase  lag  as  well  as 
amplitude  decrease.  Since  this  equation  is  linear,  a one- 
dimensional velocity  field  will  be  considered  and  the  effects 
of  self-diffusion  and  gravity  can  be  obtained  by  a separate 
analysis.  The  resulting  expression  for  the  particle  motion  is 


^up  q du  , „ , \ . o ,3a. 1/2 

dt  = B at  + ae  (u  - y + e (~ : 1 

where  a = ^ and  0 = -+~^T 

This  equation  was  solved  by  Chao  [1964]  using  Fourier  transform 
techniques,  which  allows  the  time  dependent  character  of  the 
velocity  field  to  be  expressed  as  a frequency  fluctuation  (w) . 
The  Fourier  transform  of  the  velocity  is 

oo 

u(w)  = / u ( t ) exp  (-i  cot)  dt.  (4.1.5) 


lt  If  (u  ' V at 

tr  /t-T 

(4.1.4) 


The  solution  for  the  particle  velocity  as  a function  of  the 
fluctuation  frequency  is 


{1  + [ 1 . 5u)/a]  * + 
U +[1.5u/a]i 


i(a)/a  + [1 . 5(i)/a]  ^ ) } u (i») 
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(4.1.6) 
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By  multiplying  the  velocities  u and  u by  their  respective 
complex  conjugates  and  introducing  the  energy  spectrum  density 
function  F(w),  a relation  is  obtained  between  the  intensity  of 
the  particle  motion  and  that  of  the  surrounding  fluid. 
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The  ratio  >/<u2>)  can  be  used  to  infer  the  ability  of 

the  particle  to  move  with  an  accelerating  fluid  element.  This 

ratio  is  presented  in  Figure  17  as  a function  of  the  oscillation 

frequency  w for  the  conditions:  d = 0.25,  0.5,  1,  2pm),  = 

lgm/cc,  p =*  p . (STP ) ) . One  can  infer  from  the  above  results 
ci  i r 

that  a 0.5  mm  particle  will  tollow  the  t low  with  a time  scale  of 
the  order  of  one  micro-second.  Using  the  data  of  Wyatt  at 
x/M'10,  the  Kolmogoroff  time  scale  (v/t  ) has  a minimum 
value  of  10  3sec.  Hence,  the  .5pm  particle  is  expected  to 
follow  the  most  rapid  fluctuations  to  be  found  in  the  physical 
flow  field. 

Further  insight  into  the  aerodynamic  drag  effects  can  be 
gained  by  considering  the  behavior  of  a particle  in  a swirling 
flow  which  simulates  the  motion  of  the  small  scale  eddies. 
Hirschkron  and  Echrich  11964]  have  analyzed  the  motion  of  a 
particle  in  an  irrotational  vortex  (u(1  = c^/r)  . The  salient 
features  of  this  analysis  are  presented  in  the  following. 

A spherical  particle  is  assumed  and  the  analysis  is  re- 
stricted to  the  low  Reynolds  number,  "Stokes  Flow"  regime.  The 
drag  coefficient  cD  is  given  as 

= 12 (4.1.8) 

••  • D i,-*  -* , | 

P I [u-v]  I a 


where  u is  the  fluid  velocity,  v is  the  particle  velocity,  and 
a is  the  particle  radius.  The  expression  for  the  drag  force 
F is 


? = ^ cD  p (u  - v)  • (u  - v)  na^ 


m = 6nua  [u  - v] 


(4.1.9) 


The  fluid  velocity  field  is  assumed  to  be  that  of  an  irrotational 


vortex : 


u = 0 


(4.1.10) 


The  acceleration  of  the  particle  can  be  written  in  terms  of  the 

spatial  variations  of  the  velocity  field  as 
...  ..  2 


dv 
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where  p = (6iuia/m  ) . 


These  equations  were  solved  numerically  and  the  results 

were  plotted  in  normalized  form.  The  radial  position  (r/r  ) 

2 ° 

of  the  particle  as  a function  of  time  (t  C /r  ) is  of  interest 
r-  v c 

for  the  small  scales  of  the  motion.  Choosing  an  eddy  of  the 
Kolmogoroff  size 


r(  = (v^/e)1^^  = 0.123mm 


(4.1.12) 


2 3 -5  2 

where  c = 14.7  m /sec  at  t = 0.04  sec  for  v = 1.5  x 10  m /sec. 

The  r value  is  evaluated  from  the  data  of  Wyatt  [1955] . The 

kinetic  energy  of  the  dissipative  eddies  is  approximated  by 

q = (vc)1^:  consequently,  a characteristic  value  of  u(1  is 

given  by  u()(r  )-  0.122  m/sec.  The  above  value  for  u^(r  ) 

allows  Cv  to  be  evaluated  as:  Cv  = 15  mm'Vsec.  Making  use 

of  the  Hirschkron  and  Echrich  solution  for  a 0.5-pm  particle 

in  air,  r/r{  does  not  exceed  1.2  for  an  elapsed  time  of  1.0 

second.  The  experimental  times  of  interest  are  shorter 

than  this;  it  is  again  inferred  that  the  particle  suitably 

follows  the  fluid  motion. 
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4.2.  Gravitational  Effects 

The  terminal  velocity  is  calculated  usinq  the  Stokes  draq 
relationship  givinq 

(4.2.1) 
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For 


a 0.5-um  particle  with  specific  gravity  1.0,  V^,  = 0.05  mm/sec. 
The  distance  travelled  by  the  particle  in  the  elapsed  time  of 
the  experiment  is  an  order  of  magnitude  smaller  than  the  Kolmogoroff 
scales . 

4.3.  Particle  Diffusion  by  Molecular  Agitation 

The  random  molecular  motion  in  the  neighborhood  of  the 
particle  will  cause  it  to  move  away  from  the  spatial  volume  of 
its  origin.  In  the  presence  of  a fluid  velocity  (i.e.,  a drift 
velocity  of  the  molecules)  this  motion  would  be  with  respect  to 
the  convected  spatial  volume  which  is  defined  at  the  initiation 
of  the  experiment.  The  characteristics  of  the  motion  can  be 
assessed  by  kinetic  theory  considerations  as  shown  in  the 
following . 

The  root  mean  square  velocity  of  a particle  of  mass  m is 
given  in  terms  of  the  Boltzman  constant  and  the  absolute 
temperature  T as 


1/2 


3KrT 

v = ( —5-  ) 

m. 


(4.3.1) 


For  a 0.5-vim  particle  of  specific  gravity  1.0  and  for  300°K, 
v = 4.6  mm/sec.  This  velocity  is  small  (£  1 percent)  with 
respect  to  the  rms  continuum  fluid  velocities  of  interest  in 
the  present  experiment. 

The  displacement  of  the  fluid  particle  can  be  evaluated 


from  the  diffusion  coefficient,  D 


pression 


D = 
P 


v; 
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It  is  defined  by  the  ex- 

(4.3.2) 


For  the  conditions  of  a 0.5-vim  particle  (a  = 0.25  vim)  with 
300°K  air  as  the  host  fluid, 


D = 
P 


6.7  x 10  ^ m^/sec 
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In  order  to  relate  to  the  displacement  of  the  marking  particles, 
consider  a field  of  particles  in  which  the  particle  concentration 
c is  initially  given  as 

c(z,o)  = cQ  for  z <_  o 

c(z,o)  = o for  z > o 

and  in  which  the  particle  concentration  c follows  the  Fickian 
law  of  diffusion: 


i£  . D 

3t  >>  a*2 


(4.3.3) 


The  resulting  concentration  field  for  D = constant  is  qiven  as 
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(4.3.4) 


A width  measure  for  the  migration  of  the  particles  can  be 
specified  as 


6=4  vD  t 

P P 


(4.3.5) 


For  the  conditions  used  above,  6 = 0.047  mm  for  an  elapsed 

u t P 

time  of  1.0  sec  (— ~ 250).  Hence  the  migration  of  the 
particles  in  one  second  is  0.134  times  the  scattering  volume 
dimension.  This  gives  credence  to  the  assumption  that  the 
contaminant  is  "non-dif fusive" . 

4.4.  influence  of  Particles  on  Carrier  Gas  Properties 

The  experimental  technique  of  marking  a given  spatial 

volume  requires  a substantial  number  (=20)*  of  marking  particles. 

The  influence  of  these  particles  on  the  carrier  (or  host)  gas 

is  considered  in  the  following. 

£ 

For  a concentration  of  10  , 0.5-qjm  diameter  particles/cc , 
the  density  increase  with  respect  to  the  uncontaminated  air  is 
0.1  percent.  The  influence  on  the  viscosity  can  be  assessed 
from  the  Einstein  equation;  (see  Hirschf elder , Curtiss  and 
Bird  (1954])  u ~ p(l  + 4*)  where  <p  is  the  volume  ratio.  The 


*The  particles  may  be  assumed  to  follow  a Poisson  distribution 
as  regards  their  distribution  in  the  spatial  volumes.  Hence, 
the  standard  deviation  of  the  particle  number  per  element  is 
N * where  N is  the  average  value.  Consequently,  N should  be 
large  to  maintain  an  unambiguous  scatterinq  signal. 
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change  in  viscosity  as  a result  of  the  additional  particles 

-4 

is  10  percent. 

The  influence  of  the  particles  on  the  hydrodynamic  sta- 
bility of  the  host  fluid  has  been  considered  by  Saffman  (1961] . 
His  results  can  be  used  to  infer  that  there  is  a negligible 
influence  on  the  stability  resulting  from  the  presence  of 
these  0.5-um  diameter  particles  in  air. 
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5.  RESULTS 

5.1  Flow  Visualization 


The  motion  of  the  gas  inside  the  mixing  chamber  was  photo- 
graphically recorded;  a significant  feature  of  the  initial 
motion  was  revealed  by  this  procedure.  Marking  particles, 
with  a concentration  of  106  particles/cc , were  initially  intro- 
duced in  the  upper  volume;^they  were  illuminated  by  a light 
sheet  formed  with  a cylindrical  and  spherical  lens  system. 

The  motion  pictures  were  taken  at  50  frames/sec;  the  exposure 
time  was  150  ^ sec.  The  first  visible  frame  could  not  be 
precisely  correlated  with  the  passage  of  the  grid;  it  repre- 
sents an  experimental  time  of  0-0.02  seconds.  Representative 
frames  are  shown  in  Figure  18.  The  black/white  (or  marble 
cake)  mixing  character  of  the  particle  gas  is  apparent  albeit 
some  resolution  is  undoubtedly  lost  with  the  grain  of  the  film 
and  the  relatively  long  exposure  time. 

The  first  frame  of  Figure  18  is  considered  to  be  quite 
significant.  Specifically,  the  "fuzzy"  end  of  the  marked 
particles  in  the  0-0.02  sec.  time  span  was  not  anticipated; 
it  suggests  a rapid  mixing  in  the  neighborhood  of  z = 0 which 
precedes  the  large  scale  mixing  induced  by  the  grids.  The 
source  and  character  of  this  initial  mixing  process  has  been 
determined. 

The  boundary  layers  which  grow  upon  and  which  are  shed 
from  the  splitter  plate  form  a wake  of  positive  and  negative 
vortex  sheets.  In  addition,  the  regular  vortex  structures 
shed  from  the  horizontal  cylinders  of  the  biplane  grid  will 
create  an  organized  perturbation  of  the  vortex  sheets.  The 
antisymmetric  part  of  the  perturbation  leads  to  an  instability 
of  the  vortex  sheets  with  a subsequent  roll  up  or  amalgamation 
of  the  vorticity  into  vortex  clouds.  Abernathy  and  Kronauer  [1962] 
have  computed  the  behavior  of  these  vorticities.  Their  non- 
linear analysis  shows  that  the  vortex  sheets  are  unstable  to 
any  wavelength  disturbance  and  that  the  clouds  may  be  comprised 
of  opposite  sense  vorticity  with  a consequent  reduction  in  the 
net  strength  of  the  vortex  cloud. 

+ The  particles  were  introduced  into  the  lower  half  volume  for 
the  quantitative  experiments. 
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The  spacing  of  the  regular  vortex  structures  shed  from 
the  horizontal  rods  can  be  determined  from  the  data  of 
Roshko  [1954]  and  the  conditions  of  the  present  experiment. 
The  pertinent  Reynolds  number  is 

U d 

Rn  = ~ v-  = 2000  (5.1.1) 

and,  from  Roshko,  the  corresponding  Strouhal  number  is 


fAd 

ST  = U = * °"2 


(5.1.2) 


The  wavelength  of  the  disturbance  can  be  calculated  from  f,  as 


A = U /f.  = d/0.2 

O A 

~ 2.54  cm 


(5.1.3) 


From  Schlicting  (1968),  the  impulsively  started  boundary 
layer  6 grows  as 


6 = 4/vt  (5.1.4) 

The  elapsed  time  t is  the  time  required  for  the  grid  to  move 
one  half  of  its  total  travel  across  the  box;  the  resulting 
value  for  6 is  a2.54  mm.  The  splitter  plate  is  approximately 
1.56  mm  thick  and  assuming  that  the  laminar  boundary  layer 
vorticity  is  concentrated  at  ( ~ 6/3 ) , the  h value  is  z3.26  mm. 
Hence  h/A~0.13.^ 

Abernathy  and  Kronauer  present  the  formation  of  the  vortex 
clouds  for  h/Aj>0.12.  As  shown  by  the  figure,  the  "mixing"  is 
quite  advanced  for  a unit  value  of  the  dimensionless  time 
tU'/'\#  where  U'  is  the  velocity  parameter  introduced  by 
Abernathy  and  Kronauer.  (U'=  plate  velocity/2.)  Consequently, 
the  time  to  achieve  significant  mixing  (or  dispersion  of  the 
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Significantly,  this  elapsed  time  is  quite  small  with  re- 
spect to  the  time  scales  which  characterize  the  subsequent 
mixing;  specifically,  the  minimum  Kolmogorov  time  in  the 
homogeneous  turbulence  field  is  * 5 msec.  Also,  the  7 msec  to 
establish  the  initial  mixing  layer  is  comparable  with  the 
elapsed  time  (3.5  msec)  between  two  successive  scans.  Conse- 
quently, the  experimental  technique  will  adequately  identify 
the  structure  of  this  mixed  region  once  it  is  formed,  but  it 
will  not  be  able  to  document  its  evolution. 

The  analysis  of  Abernathy  and  Kronauer  show  that  the  mixed 
region  is  four  to  five  times  the  original  separation  distance 
h.  Consequently,  the  initial  conditions  for  the  quantitative 
evaluation  of  the  mixing--as  controlled  by  the  large-scale 
motions--will  be  a lower  region  of  contaminated  gas  and  an 
upper,  uncontaminated  gas,  region  separated  by  a well  stirred 
region  of  height  ~0.5M. 

The  smaller  kinematic  viscosity  of  the  Freon  gas  suggests 
that  the  initial  displacement  of  the  two  vortex  sheets  will 
be  smaller  than  the  displacement  for  the  air  case.  The  impul- 
sively started  boundary  layers  will  be  reduced  by  the  multi- 
plicative factor  [v_/v  . ] 1//2  ^ ..  ^ .. 

F air  ~27T.  This  suggests  that  the 

initial  width  measure  for  the  Freon  results  should  be  smaller; 
this  point  is  addressed  in  the  next  section. 

1 

The  later  times  in  Figure  20  and  21  show  the  mixing  for 
longer  time  periods.  It  is  evident  from  these  pictures  that 
the  interface  is  highly  convoluted  and  stretched  by  the  turbu- 
lent strain  field. 
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5.2  Data  Base  for  the  Particle  Experiments 

The  "non-dif fusive  mixing"  experiments  are  achieved  by 
contaminating  the  lower  volume  (gas  B)  with  a cloud  of  submi- 
cron particles  whose  dimensions  are:  (1)  sufficiently  large 
that  they  are  not  significantly  influenced  by  the  Brownian 
motion  of  the  surrounding  gas  molecules  and  (2)  sufficiently 
small  that  their  "slip  velocity"  with  respect  to  the  continuum 
motion  is  negligible.  The  nominal  diameter  of  the  particles, 
viz.,  0.5ym  satisfies  these  joint  criteria.  The  turbulent 
transport  in  the  chamber  is  revealed  by  the  presence  of  B in 
the  upper  half  volume  and  that  of  A in  the  lower  volume.  Par- 
anthetically , it  can  be  noted  that  the  incompressible  nature 
of  the  mixing  experiment  requires  that  the  volumes  of  the  two 
gases  are  invarient  with  respect  to  time;  however,  this  fact 
cannot  be  readily  used  in  the  interpretation  of  the  data  since 
information  from  only  one  spatial  line  is  available  and  because 
sub-scattering  volume  mixing  may  be  present  which  would  yield 
an  excessive  estimate  for  the  presence  of  B.  In  principal, 
the  presence  of  B in  the  scattering  volume  is  easily  detected. 
The  light  scattering  from  the  particles  (of  B)  creates  a photo- 
multiplier tube  output  signal  which  is  many  times  that  for  the 
noncontaminated  gas.  A comparator  circuit  readily  distinguishes 
the  two  states  (a  nominal  signal/noise  ratio  _>  10  was  achieved) 
and  creates  a binary  pulse  train  with  a "0"  indicating  the 
presence  of  A and  a "1"  indicating  the  presence  of  B. 

The  electronics  processing  circuitry  converted  the  pulse 
train,  from  the  comparator  circuit,  into  data  words  which  were 
compatible  with  the  TI  minicomputer.  The  data  was  then  trans- 
ferred to  digital  magnetic  tape  for  further  processing  on  the 
CDC  6500  central  computer. 

The  experiments  were  carried  out  with  air  as  the  working 
medium  for  69  realizations  and  with  Freon  12  for  100  realiza- 
tions. The  marking  particles  were  placed  in  the  lower  half 
volume  for  all  cases.  A 74  mm  scan  was  provided  by  the  mimi- 
fication  of  the  optical  system  and  there  were  205  readings  per 


41 


scan;  hence,  each  data  point  was  sampled  from  some  fraction  of 
a 0.36  mm  section  of  the  laser  beam  in  the  mixing  chamber.  The 
nominal  scattering  volume  height  (h)  may  be  estimated  from  the 
slit  width  (w) , the  magnification  ratio  (M)  and  the  angle  of 
intersection  between  the  slit  and  the  vertical  laser  beam  (<J»)  . 
Specifically,  h = lw/M  cos  = [0.25/1.36  cos(tt/4)]  = 0.26  mm. 
The  actual  probe  volumes  were  both  somewhat  larger  and 
smaller  than  this  nominal  value;  however,  less  than  1 percent 
of  the  scattering  volumes  exceeded  the  0.36  mm  spacing  noted 
above  and  consequently  they  can  be  considered  to  be  non- 
overlapping. The  scans  were  executed  in  3.11  ms;  the  period 
of  the  scans  was  3.5  ms.  The  designated  scan  number  is 
precisely  known  from  the  data  record;  however,  its  time, 
relative  to  the  passage  of  the  grid,  is  uncertain  since 
there  was  no  coupling  between  the  initiation  of  the  scan  and 

the  grid  motion.  Consequently,  the  designated  time  should 

U t 

be  understood  as  = 0.433. 

The  data  base  which  results  from  these  experiments  is  re- 
presented by  the  $/blank  map  shown  in  Figure  15.  The  presence 
of  gas  B,  as  marked  by  the  particulate  matter,  is  represented 
by  the  presence  of  the  unmarked  fluid  A,  is  signified  by 

a "blank"  space.  The  map  represents  data  for  every  other 
location  (up  to  location  200)  in  keeping  with  the  lateral 
dimension  from  the  line  printer.  The  quantitative  results  make 
use  of  the  complete  set  of  205  locations.  The  interpretation 
of  the  experiments  requires  a careful  examination  of  numerous 
aspects  of  the  physical  processes  in  the  chamber  and  the  details 
of  the  data  acquisition  process;  these  matters  are  dealt  with 
in  the  following  sections.  Similar  maps  for  stable  and  un- 
stable cases  involving  Freon  12  and  air  are  presented  in 
Appendix  C.  A map  showing  a complete  set  of  240  scans  for  an 
air/air  case  is  also  presented  in  this  appendix. 

5.3  Ensemble  Average  Distributions 

The  experimental  results  for  the  ensembles  of  69  and  100 
realizations  of  the  air  and  the  Freon  12  gases  are  considered 
in  this  section.  The  instantaneous  concentration  of  gas  A or 
B at  a given  location  and  at  a given  time  from  the  initiation 
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of  the  experiment  is  termed  T(z,t);  for  the  non-dif f usive  ex- 
periments its  possible  values  are  0 or  1.  The  ensemble  average 
values  of  T are  identified  as  <T(z,t)>  and  such  values  have 
been  evaluated  for  a discrete  set  of  times  corresponding  to 
scans  (1,  2,  4,  6,  8,  11,  14,  17,  21,  25,  29,  33,  37,  41,  60, 

80,  100,  120,  140,  160,  180,  200,  220). 

Since  the  mixing  field  for  the  neutral  density  condition 
is  at  least  expected  to  be  symmetric,  and  since  the  turbulent 
transport  effects  in  a homogeneous  turbulence  field  are  known 
to  be  readily  approximated  as  a Gaussian  process,  a normal 
distribution  was  examined  as  a way  to  characterize  the  lateral 
distribution  of  the  ensemble  average  values.  Analytical  con- 
siderations regarding  the  use  of  the  normal  distribution  have 
been  presented  in  Section  3.2.  The  data  were  found  to  follow 
such  a distribution,  at  least  for  a first  approximation,  and 
the  subsequent  evaluation  of  the  data  will  rely  upon  the  para- 
meters of  the  normal  distribution  to  provide  reference  quantities 
for  the  spatial  distributions  of  the  average  concentrations. 

A representative  data  set  is  presented  in  Figure  22. 

The  dimensional  plot  <T(z,t)>  is  shown  in  part  (a);  the  solid 
curve  of  this  figure  is  the  normal  distribution  found  by  a 
linear  least  squares  fit  to  the  representation  of  the  data  on 
the  "probability  coordinates"  of  part  (c) . The  probability 
coordinates  are  "defined"  in  part  (b) . The  standardized  normal 
cumulative  distribution  function,  defined  as 

~2 

A(£)  = exp  (-  \)  da  (5.3.1) 

' ' —CO 

see  Dubes  [1968],  is  shown  on  part  b.  The  coordinate 
transformation 

5 * % ~l  <*  - *c>  (5-3.2) 


related  the  standard  Gaussian  distribution  to  the  ensemble 
average  concentration  distribution,  <T(z,t)>.  The  symbol  zc 


A3 

represents  the  centerline  of  this  distribution;  i.e., 

<F(z  ,t)>  =0.5.  o represents  the  standard  deviation  of  the: 
c m 

spatial  distribution;  i.e.,<F(z  * o ,t)>=  0.5*0.34.  The  stan*- 

c m 

dard  deviation  represents  a convenient  width  measure  of 

the  distribution.  Figure  23  represents  the  ensemble  average 

data  <F(z,t)>  and  Figure  24  shows  the  resultant  transformation 

onto  probability  coordinates  for  Freon  and  air  at  UQt/M  = 12. 

Note,  that  zc  is  different  from  the  centerline  of  the  mixing 

chamber  (z^) . Apparently,  this  is  a result  of  large  scale 

gross  transport  of  the  mixing  region.  It  is  assumed  that  this 

transport  does  not  directly  effect  the  turbulent  mixing  near 

the  centerline  from  which  o is  calculated.  The  essential 

m 

features  of  the  ensemble  average  data  are  represented  by 
straight  lines  on  probability  coordinates;  see  Figures  25 
and  26  for  a composite  representation  of  the  air  and  Freon  12 
results.  The  standard  deviation  between  the  cumulative  dis- 
tribution function  and  the  experimental  data  are  also  shown 
in  the  tables  included  with  these  plots. 

The  width  measures,  o , for  the  Freon  and  air  data  are 

m 

shown  in  Figure  27.  This  figure  also  shows  the  physical 
size  of  the  round  rod  and  the  mesh  size  of  the  biplane  grid. 
The  evolution  of  the  mixing  field  within  the  chamber  is 
shown  in  Figure  28.  The  spatial  locations  of  <F>  = 0.15  and 
<F>  = 0.85  are  used  to  mark  the  "boundaries"  of  the  mixing 
field  as  a function  of  the  elapsed  time  from  the  initiation 
of  the  experiment  (UQt/M). 

5.4  Measures  of  the  Instantaneous  Scans 

The  individual  measurements  of  F(z,t)  in  the  "instanta- 
neous" scans  of  the  turbulent  mixing  field  represent  a very 
considerable  amount  of  information  which  is  only  partially 
characterized  by  the  ensemble  <F(z,t)>  distributions. 

Insight  into  the  detailed  structure  and  processes  of  the 
mixing  field  can  be  gained  if  suitable  measures  of  the 
instantaneous  scans  can  be  formulated.  One  such  formulation 
has  been  achieved;  it  is  described  below. 
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Consider  two  extreme  cases  for  the  instantaneous  concen- 
tration field;  both  are  capable  of  yieldinq  the  observed  <T(zrt)> 
following  the  ensemble  average  process.  The  two  cases  are 
introduced  in  the  form  of  statements.  "The  instantaneous 
concentration  field  is  such  that  an  instantaneous  scan  re- 
veals : 

Case  1.  a well  mixed  field  with  transverse  dimension  and  conr- 
centration  r(z,t)  "similar"  to  that  of  vr(z,t)N." 

Case  2.  a relatively  thin  mixed  region  (i.e.,  one  which  is 
small,  for  example,  with  respect  to  a (t))  which  is 
laterally  transported  by  the  convective  action  of  the 
large  scale  turbulent  motion  to  yield  the  observed 
<r(z,t)N  distribution." 

It  is  not  expected  that  either  of  these  models  will  be 
supported  in  their  extreme  condition;  however,  they  do  serve 
as  a guide  for  the  definition  of  measures  to  evaluate  the 
character  of  the  instantaneous  scans.  Specifically,  defini- 
tions for  i)  the  center  of  the  instantaneous  scan  and  ii)  the 
width  of  the  mixed  region  about  the  instantaneous  center  are 
required.  The  two  definitions  are  separately  considered  in 
the  following;  the  center  of  the  instantaneous  scan  is  defined 
by  the  following  operations. 

a. )  Determine  the  number  of  blanks  (or  absences  - A) 

which  signify  the  presence  of  fluid  A in  the 
scattering  volumes  and  designate  this  number  as  NA. 
Similarly,  identify  the  number  of  filled  regions 
and  designate  this  number  as  NB  (NB=205-NA) . 

b.  ) The  center  of  the  instantaneous  scan  will  be 

designated  as  the  position  along  the  scan  where 
the  number  of  scattering  volumes  between  the  upper- 
most position  of  the  scan  and  the  instantaneous 
center  is  equal  to  NA. 

Note  that  this  definition  implicitly  assumes  that  the  B 
fluid  does  not  extend  beyond  the  upper  limit  of  the  scan  and 


that  the  A fluid  does  not  extend  beyond  its  lower  limit.  Also 
note  that  the  upper  and  lower  designations  are  understood  in 
terms  of  the  physical  experiment;  in  terms  of  the  $(E  B)  map 
of  Figure  15,  "uppermost"  corresponds  to  "leftmost".  A 
schematic  representation  of  the  above  operations  is  provided 
in  Figure  29. 

The  position  of  the  instantaneous  center  for  the  ith  scan 

and  for  jth  experiment  will  be  designated  as  z . ) ..  The  lo- 

cation  of  the  average  of  the  instantaneous  centerlines  ^z  . > 

ci 

and  its  displacement  from  the  ensemble  average  centerline, 

<zci>  - <zc>  (see  Figure  29),  are  of  interest.  The  statistical 

characteristics  of  the  z .)  . - < z . > values  allow  one  feature 

ci  j ci 

of  the  two  cases  to  be  assessed;  namely,  if  the  variation  of 
the  centerline  position  is  large  (with  respect  to  a for 
example)  then  support  for  case  2 is  inferred.  The  standard 
deviation  of  the  z^)^  - <zc^>  population  is  defined  as 
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(5.4.1) 


where  N is  the  total  number  of  experiments.  The  quantity  o^z 
is  compared  with  om  in  Figure  30a  and  30b;  the  interpreta- 
tion of  these  results  is  discussed  in  the  next  section. 

The  width  of  the  mixed  region  about  the  instantaneous 
scan  can  be  assessed  in  a manner  quite  similar  to  that  of  the 
width  measure  for  the  ensemble  averages.  (This  will  offer  the 
subsequent  advantage  of  a basis  for  the  comparison  of  their 
relative  magnitudes.)  Specifically,  the  emsemble  average  of 

the  T field  with  respect  to  z . ) was  formed  as 

ci  j 

<r(z  - zci)>  = £ E T(z  - Z .)  ) (5.4.2) 

1 j = l J 

where  n<N  and  n represents  the  number  of  experiments  which 
possess  a measured  value  of  I'  at  the  z location  in  question. 
The  <f (z  - zc ^ ) > so  evaluated  were  then  fit  with  the  normal 


distribution  as  described  by  equations  5.3.1  and  5.3.2.  The 
resulting  value  of  the  width  measure  is  termed  the 

values  of  a ^ i for  the  69 (air/air)  and  the  100 (Freon  12/Freon 
experiments  are  compared  with  the  width  measure  of  the  en- 
semble average  scans,  o )i#  in  Figures  30a  and  30b.  These 
results  are  discussed  in  the  next  section. 
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6.  DISCUSSION  OF  RESULTS 

6.1.  The  Ensemble  Average  Concentration  Distributions 

Figures  25,  26,  and  27  represent  the  primary  documenta- 
tion of  the  experimental  results  for  the  ensemble  average 
concentration  distributions.  These  figures  include  a tabula- 
tion of  the  standard  deviation  (S.D.)  values  which  indicate 
the  validity  of  representing  the  ensemble  average  data  with 
a Gaussian  Cumulative  Distribution  Function  (C.D.F.).t  The 
S.D.  values  are  computed  from  the  relationship  between  <T(z,t)> 
and  the  C.D.F.  in  the  range  0.03  £<T(z,t)>£  0.98.  These  limits 
reflect  the  relatively  small  sample  size,  and  hence  the  limited 
resolution,  of  the  experimental  data  (1:69  air  and  1:100 

Freon  12)  and  are  selected  such  that  the  o and  z values  are 

m c 

"minimally"  effected  by  the  selection  of  an  arbitrary  segment 

of  the  experimental  data.  (Trial  runs  on  one  data  set  show 

only  a nominal  influence  of  the  cut-off  values  on  o ; speci- 

m 

fically,  c>m  for  UQt/M  = 18  did  not  change  by  more  than  2fc 
when  cut-offs  representing  0.014  - 0.986,  0.028  - 0.972,  and 
0.043  - 0.957  were  examined.)  The  relatively  small  values 
of  S.D.  (compared  to  unity)  suggest  that  the  Gaussian  C.D.F. 
provides,  at  least,  a satisfactory  representation  of  the 
experimental  data.  Since  this  analytical  form  is  also  sup- 
ported by:  (i)  the  (assumed)  Gaussian  displacement  of  fluid 
elements,  see  [3.2.6],  and  (ii)  (independently)  by  an  assumed 
"Dt  = constant"  model  of  the  turbulent  transport,  see  13.2.5], 
this  compact  presentation  of  the  <I'(z,t)>  results  is  both 
efficient  and  apparently  justified. 

The  experimental  data  reveal  that  the  values  for  the 
air  experiments  are  larger  than  the  o values  for  the  Freon  12 
experiments  at  a given  UQt/M.  This  result  was  not  expected; 

tThe  experimental  data,  <l'(z,t)>,  are  used  to  define  the  Gaussian 
distribution  by  the  selection  of  z and  o which  minimizes 
(least  squares)  the  difference  between  the  experimental  data 
values  and  the  C.D.F.  values.  The  latter  is  defined  as 

(I?)4o  / exp[  - J [-j-S]  da. 

m -<»  m 
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it  is  considered  to  be  quite  significant.  If  the  large  scale, 

energetic  motions  in  the  turbulence  field  are  assumed  to  be 

independent  of  the  Reynolds  number,  then  the  smaller  a values 

m 

for  the  Freon  12  experiments  suggests  that  the  transport  is 
not  solely  controlled  by  the  large  scale  motions.  Additional 
interpretations  of  this  result  are  made  possible  by  the  ana- 
lytical considerations  of  Articles  6.2  and  6.3. 

The  spatial  evolution  of  the  mixing  region,  which  is 
suggested  by  these  results,  is  presented  in  Figure  28.  it 
is  of  interest  to  note  that  the  boundaries  of  the  ensemble 
average  mixing  region  are  different  for  the  air  and  Freon  12 
experiments  and  that  gross  differences  are  observed  for 
UQt/M  > 50.  It  is  inferred  that  subtle  differences  in  the 
two  sets  of  experiments  are  probably  responsible  for  the 
differences  in  the  spatial  evolution.  It  is  also  inferred, 

from  the  relatively  self-consistent  behavior  of  o , oT  and 

ml 

°Az  ' t^iat  fcViese  differences  are  not  significant  in  the  in- 
terpretation of  the  experimental  results.  The  Freon  12  results 

show  a distinct  trend,  beginning  at  approximately  Uot  _ 7q_q5 

M ' 


where 


the  entire  mixing  field  appears  to  be  transported  to  small  z 
values.  Since  the  extent  of  the  field  that  is  available  to  be 
fit  by  the  Gaussian  C.D.F.  is  reduced  below  a marginal  value,  the 
om  for  UQt/M>110  are  not  presented  for  the  Freon  12  carrier  gas. 

An  interesting  comparison  can  be  made  between  the  tur- 
bulent diffusivity  represented  by  and  the  molecular  dif- 

fusivity  D of  the  carrier  gas.  Equation  (3.2.9b)  relates 

2 

the  turbulent  diffusivity  io  the  slope  of  the  z^  (t)  curve. 

The  following  results  are  obtained. 

2 

Dt (air)  = 14-6  cm  /sec 

for  5£  Uot/M  £ 20  (6.1.1a) 

Dt (Freon)  = 12.6  cm^/sec 


Dx(air  and  = 29.2  cm^/sec  for  20<  Uot/M  <200  (6.1.1b) 

Freon) 


The  following  ratios  are  obtained  for  the  turbulent  dif- 
fusivity to  the  molecular  diffusivity  for  air  in  this  experi- 


ment. 
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(6.1.2a) 


73  5 < UQt/M  <20 

t 

146  20<  U t/M  <200  (6.1.2b) 

6.2.  Measures  of  the  Instantaneous  Scans 

Figures  30a  and  30b,  which  present  the  results  from  the 
instantaneous  scan  evaluations  for  air  and  Freon  respectively, 
strongly  support  the  "well  mixed"  condition  described  by 
Case  It  of  Article  5.4.  Specifically,  the  values  are  with- 
in 10  percent  of  the  o values  at  a given  U t/M.  The  o. 

m ^ o Az 

values  for  air  and  Freon,  which  provide  a measure  of  the 
range  covered  by  the  instantaneous  centerlines,  are  substan- 
tially smaller  than  those  for  a ^ (or  o ) and  they  show  a much 
slower  growth  following  their  rapidly  established  "initial"  value. 

Only  the  instantaneous  scans  for  U t/M  <_  35  are  considered; 
the  ensemble  average  results  suggest  that  within  this  bound, 
the  mixing  field  is  confined  within  the  74  mm  scan  length. 

There  is  clearly  more  "scatter"  in  the  air  results  than  in 
the  Freon  results.  Presumably,  this  is  related  to  the  number 
of  experiments  which  comprise  the  ensemble:  air:  69,  Freon:  100. 
The  Gaussian  C.D.F.  provides  a satisfactory  approximation  to 
the  <T(z  - zc^) > data  for  both  the  Freon  and  the  air  experi- 
ments. This  inferrence  is  based  upon  the  acceptably  small 
values  of  the  standard  deviation  (S.D.)  between  the  <T(z  - zci)> 
experimental  results  and  the  <T(z  - zc^)>  values  computed  from 
the  appropriate  C.D.F.  (The  latter  is  characterized  by  the 
Oj  value  described  in  Article  5.4.  and  is  determined  from  a 
least  squares  fit  to  the  data.  Note  that  the  S.D.  is  based 
upon  the  population  [<T(z  - 2C^)>  - <T(z  - z ^)>  ] where 

data,k  Gaus.  C.D.F. ,k 

k is  the  index  number  representing  the  discrete  sample  volume. 

1 < k < 205. 

tCase  1:  "The  instantaneous  concentration  field  is  such  that 
an  instantaneous  scan  reveals  a well  mixed  field  with  trans- 
verse dimension  and  concentration  T(z,t)  similar  to  that  of 
<T(z,t) >. " 


D (air) 
Dt 

DTair) 
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The  o^z  values  for  the  Freon  12  and  air  experiments  are 
essentially  identical;  this  comparison  is  shown  on  Figure  30b,. 
This  is  considered  to  be  a most  significant  result.  It  was  ! 
expected,  a priori,  that  the  experimental  results  which  were 
dependent  upon  the  large  scale  motions  would  not  exhibit  a 
dependence  on  the  Reynolds  number.  The  o^z  values  are  con- 
sidered to  singularly  represent  the  large  scale  transport  of 
the  interfacial  region;  hence,  they  are  expected  to  be  the 
same  for  the  Freon  12  and  air  experiments.  The  agreement 
between  the  two  set  of  experiments  is  considered  to  be  a con- 
firmation of  the  experimental  technique  and  the  above  hypo- 
thesis. However,  it  could  (naively)  be  assumed  that  o is 

m 

also  singularly  dependent  upon  the  large  scale  motions.  This 
(plausible)  assumption  is  in  conflict  with  the  observations 
for  om  and,  as  shown  in  Figure  30,  for  Oj  as  well.  Specifi- 
cally, oI(air)  > OjfFreon  12).  Theoretical  considerations, 
which  are  based  upon  the  Lagrangian  dispersion  analyses  of 
Appendicies  D and  E,  are  presented  in  Article  6.3.  It  is 
shown  in  that  article  that  a rational  explanation  of  the 
smaller  am  values  can  be  formulated.  It  is  presumed  that  the 
physical  influence  of  the  kinematic  viscosity,  which  was  re- 
sponsible for  the  relatively  lower  om  values  for  the  Freon  12 
experiments,  is  also  responsible  for  the  lower  Freon  12  Oj 
values.  Consequently,  a consideration  of  the  o^.  results  will 
be  deferred  until  Article  6.3. 

6.3  Inference  of  the  Lagrangian  Scales  in  the  Present  Experiment 

As  discussed  in  Section  3.2.,  the  dispersion  parameter 
z is  considered  to  be  numerically  equal  to  the  o value. 

For  this  experiment  the  contaminant  which  occupies  the  lower 
half  of  the  chamber  has  a physical  "release  time"  tr  which 
corresponds  to  the  instant  the  grid  and  splitter  plate  pass 
the  observation  plane.  The  theoretical  developments  of 
Appendix  D,  specifically,  the  formulation  of  the  dispersion 
using  two  separate  time  scales,  can  be  examined  in  the  con- 
text of  the  o experimental  values.  The  micro  and  macro  time 


I 


fl 


I 

L 
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scales  can  be  estimated  by  the  procedure  outlined  in  Appendix  E. 
These  applications  of  the  theoretical  results  are  carried  out; 
in  this  section. 

I 

6.3.1.  A Priori  Estimation  of  the  Lagrangian  Microscale 

Utilizing  Shlien's  experimental  data  and  the  theoretical 
development  by  Hinze  [1975],  an  estimate  of  the  normalized 
microscale,  t*,  (tl/C^)  was  obtained  in  terms  of  the  grid 

Reynolds  number  and  the  parameters  in  the  kinetic  energy  decay 
law,  see  equation  (E.7.8).  For  the  present  experiment,  these 
parameters  have  the  following  numerical  values  for  decaying 
grid  turbulence: 


AIR 


FRF.ON  12 


A = 31 

U_  = 6.28  m/sec 
o 

M = 2.54  cm 
-4  2 

v_ . = 0.15x10  m /sec 

all 

n = 1.27 


v 


A = 27 

U = 6.28  m/sec 

M = 2.54  cm 

* = 0.17  v 

freon  air 

n = 1.27 


From  Equation  (E.7.8.)  the  normalized  microscale  can  be  de- 
fined as 

t*  = (i  - 1>  in  u ♦ -V3/4  ijH'2Y  (6-3.1.) 

CL  2 m 

This  leads  to  the  following  values: 


T*  (air)  * *26 
T (Freon)  ~ 


(6.3.2.) 


2 

The  theoretical  result  for  z , presented  in  (3.2.22)  and 

P 2 

plotted  in  Figure  33,  shows  that  z^  increases  for  an  increas- 
ing value  of  x*.  Hence,  the  difference  in  the  predicted 
t*  values  for  air  and  Freon  12  suggests  that  the  dispersion, 
using  air  as  a carrier  gas,  will  be  greater  than  that  for 

the  Freon  12  carrier  gas.  The  o values  of  Figure  27  are  in 

m 

= .17  obtained  from  Shlien’s  T1972]  data  for  dispersion 
in  decaying  homogeneous  turbulence. 


qualitative  agreement  with  this  prediction;  quantitative  con-, 

i . ! 

siderations  are  presented  below. 

j 

For  the  present  experiment,  the  mixing  begins  before  the 
turbulence  field  becomes  homogeneous.  That  is,  there  is  a 
period,  UQt/M  < 5,  of  inhomogeneity  during  which  significant 
mean  velocity  gradients  and  spatially  non-uniform  velocity 
fluctuation  intensities  exist.  Sato  and  Yamamoto  [1974]  have 
quantitatively  investigated  this  period  in  a wind  tunnel  in- 
vestigation of  the  near  field  behind  a grid.  The  experimentally 
observed  results  of  the  dispersion  were  influenced  by  the 
inhomogeneous  turbulence  field  for  this  initial  time  period  and 
since  these  physical  effects  cannot  be  accounted  for  by  the 
theoretical  considerations  of  Appendix  D,  an  alternative  inter- 
pretation of  the  data  is  required. 

6.3.2.  Evaluation  of  Parameters  Which  Describe  the  Turbulent 
Dispersion 

*2 

The  analytical  structure  for  the  dispersion  z which 
was  developed  in  Appendix  D,  and  extended  in  Appendix  E,  in- 
volves the  parameters  I 

U t At . 

or  i 


In  principal,  the  available  experimental  data  can  be  inter- 
preted in  this  theoretical  context  to  infer  these  parameters. 
This  inferrential  process  must  acknowledge  that  the  analytical 
structure  is  quite  complex  and  that  the  assumptions  involved 
in  its  construction  may  invalidate  part  or  all  of  it.  It 

was  found  useful  to  create  a set  of  working  drawings.  These 

1 2 

include:  (i)  the  mean  square  particle  displacement  ^z  as 

•J*  ^ 

a function  of  UQt/M)  data  for  both  air  and  Freon  12  over 

the  full  range  of  U t/M  values  and  (ii)  an  expanded  scale 
12  ° 

plot  of  yz  for  U t/M  < 35;  see  Figures  31  and  32.  Note  that 
these  working  plot  identify  the  dispersion  as  a function  of 
the  (normalized)  laboratory  time.  The  theoretical  curves  for 


2 2 

tFrom  subsection  3.2,  z is  identical  to  o for  the  assumptions 
outlined  in  that  subsection. 


1 *2  1 *2 

jz  ( At/Atr ; x*)  and  jZ  MAt/Atr;  A*)  were  also  prepared;  see 
Figures  33  and  34.  Note  that  n = 1.27  for  the  latter  two 
plots;  this  value  is  taken  from  the  study  by  Comte-Bellot  and 
Corrsin  [1966] . 

A detailed  examination  of  these  results  reveals  that  the 
present  analytical  formulation  is  incompatible  with  the  ex- 
perimental data.  Specifically,  the  experimental  results  re- 
veal a "plateau"  in  the  z^  distribution  (25  < UQt/M  < 40)  and 
an  examination  of  Figures  33  and  34  show  that  a similar  juncture 
between  the  microscale  and  macroscale  dominated  distributions  is 

not  possible;  that  is,  the  present  analytical  formulation  only 

2 

provides  for  a juncture  of  the  two  linearly  increasing  z^ 
distributions.  Apparently,  a more  accurate  representation  of 
the  Lagrangian  correlation  function  would  be  required  to  pro- 
vide the  "plateau"  observed  in  the  experimental  data  (Figures 
31  and  32 ) . 

6. 3. 2.1.  A Posteriori  Evaluation  of  the  Microscale  and  At 

r 

The  experimental  data  can  be  used  to  evaluate  the  Lagran- 
gian micro  and  macroscales;  the  procedure  is  to  find  the 
parameter  values  which  provide  optimal  agreement  between  the 
theoretical  and  experimental  results.  This  task  is  made 
ambiguous  by  the  unknown  value  of  At^  to  use  in  the  analysis. 
Physically,  t = 0;  however,  the  analysis  is  not  valid  for 
small  times  because  of  the  strong  inhomogeneities  of  the 

turbulence  field.  Consequently,  a scheme  to  select  self- 

2 

consistent  parameter  pairs  (r*,  Atr) , which  also  yield  z^ 

values  in  agreement  with  the  data,  was  adopted.  This  scheme 

is  described  below. 

*2 

The  z (At/At  ) theoretical  results  from  the  microscale 
r *2 

analysis  (Figure  33)  show  that  z increases  linearly  for 

sufficiently  large  At/At  values.  (The  curves  are  necessarily 

+r 

terminated  at  At^/At  ) . Since  the  experimental  results  for 


+An  alternative  data  processing  scheme  which  would  remove  this 
ambiguity  is  under  consideration. 
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also  reveal  a region  of  linear  growth  (see  Figure  31) , a 
scheme  which  equated  the  theoretical  slope  to  the  processed 
data  slope  was  adopted. 

2 

The  experimental  data  for  z^  are  made  non-dimensional  by 

the  specification  of  an  arbitrary  At  value  (see  3.2.21)  and 

*2  r 

the  slope:  0.5d(z  ~ ) /d ( At/ Atr ) , is  then  determined.  A unique 
value  of  x*  can  then  be  found  which  provides  the  best  agree- 
ment between  the  experimentally  and  theoretically  determined 
slopes.  (Note  3.2.22  for  z*  (At/At  ; T*,  n) ; also  note  that 
n is  adopted  from  wind  tunnel  data  for  turbulence  kinetic  energy 

decay,  Comte-Bellot  and  Corrsin  [1966].)  This  process  yields 

2 

the  following  results  for  the  z ) data. 

p dir 

Table  II  Evaluation  of  Dispersion  Parameters 


(U0/M) 

(Atr)  : 

1.73 

2.14 

3.47 

5.17 

7.36 

10.8 

T * : 

0.23 

0.25 

0.3 

0.35 

0.4 

0.5 

(U  /M) 

(At  ) : 

3.25 

4.24 

7.91 

13.5 

22.08 

42.55 

The  Atj  values  of  this  table  are  instructive;  they  provide  a 
means  of  establishing  an  upper  bound  for  At^..  The  experimental 
data  reveal  that  the  plateau  begins  at  approximately  Uot/M  = 20. 
Hence  U AtT/M  must  be  less  than  ~ 15  since  U t /M  ~ 5.  (See 
Comte-Bellot  and  Corrsin  (1966].)  It  is  apparently  not  pos- 
sible to  select  a lower  limit  for  the  possible  combination  of 

2 

Atr  and  r*  values  since  the  z^  data  could  continue  to  grow 
linearly  beyond  the  region  which  is  dominated  by  the  micro- 
scale. However,  some  guidance  regarding  the  minimum  value 

*2 

to  be  expected  is  gained  from  the  z distributions  of  Fiqure 
32.  The  value  t*  = 0.25  has  been  adopted  as  the  plausible 
minimum.  Note  that  the  implied  juncture  region  for  this  i * 
value  would  be  quite  extensive  since  the  macroscale  dominated 
region  is  not  initiated  before  UQt/M  ~ 30. 


+From  the  expression:  At  /At  = exp [ t */ ( l-n/2 ) ] , it  is  seen 
that  At  represents  the  maximum  possible  value  of  At  for  whicfc, 
the  microscalc-dominated  analysis  could  be  used  to  describe  z 
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The  estimate  of  x * = 0.26  presented  in  6.3.1.  is,  there- 
fore, suspect.  It  is  not  clear  whether  this  prediction  would 
be  expected  to  fail  because  of  an  improper  numerical  value 
for  C or  the  inappropriate  use  of  information  from  Shlien's 

L + 

experiment.  It  is  pertinent  to  note  that  if  the  i * values 
of  (6.3.2.)  are  accepted  as  valid,  it  would  be  extremely  dif- 
ficult to  confirm  them  with  experimental  data  since  the  At^ 
values  would  be  quite  small. 

4.  4- 

The  present  experimental  data  will  be  used  to  reevaluate 
the  constant  CT  in  the  expression  (6.3.1.)  for  i*  as  a func- 
tion  of  the  parameter  A and  the  mesh  Reynolds  number  (R^) . 

For  this,  it  will  be  assumed  that  t*.  = 0.35  for  the  air 

results;  note  that  this  value  represents  the  "largest  feasible" 

2 

Atr  which  can  be  used  to  fit  the  slope  of  the  data  for 
U t/M  < 20.  For  A = 31,  R„  = 104,  and  t*  = 0.35,  CT  is 
evaluated  to  be  0.2  and  (6.3.1)  can  be  rewritten  as 

0.25 


x*  = (1  - ^)  In  (1  + 6.81  §- 

2 *M 


) 


(6.3.3.) 


This  equation  can  be  used  to  assess  the  influence  of  the  ki- 

2 

nematic  viscosity  on  the  observed  slope  of  the  z^  distribution 
between  the  air  and  Freon  12  results. 

The  calculation  of  Tprcon  recIuires  a knowledge  of  Atr. 

, + + 

It  is  reasonable  to  assume  that  At^  is  increased  by  (Uot/M  = 1) 
to  UQtr/M  = 6.17  since  t (apparent  temporal  origin  of  the 
decaying  turbulence)  for  the  larger  Reynolds  number  is  known 

to  decrease  by  approximately  this  amount;  see  Wyatt  [1955). 

• 2 
The  experimentally  observed  slope  of  the  z^  distribution  can 

tMolecular  diffusion  effects  from  the  temperature  marker  in 
Shlien’s  experiment  could  provide  a low  estimate  of  t*.  This 
is  discussed  in  Appendix  E.3. 

t+Shlien's  data  was  used  in  E.7.8.  to  evaluate  ; a value  of 

0.17  was  found.  The  present  is  approximately  18  precent 

larger  than  that  inferred  from  the  Shlien  data.  cr  ( 3/2 ) * 7 5 ; 
its  value  was  determined  to  be  4.44.  The  present  C (=  0.2) 
differs  from  CL  found  using  Shlien's  data  (i.e.,  0.17)  by 
+18  percent. 


tttThe  final  result  is  only  slightly  affected  (5%  increase)  by 
this  assumption. 


be  made  dimensionless  by  the  procedures  noted  above;  its  value 

is  0.34.  This  corresponds  to  a x*  value  of  0.24.  Hence 

Freon  12 


Freon  12 

T * 


0.24 

0.35 


= 0.69 


The  ratio  of  the  x*  values  computed  from  (6.3.3.)  using  the 
air  and  Freon  12  values  is  0.74.  The  agreement  between  these 
two  results  is  interpreted  as  a strong  endorsement  for  the 
theoretically  inferred  effect  of  the  kinematic  viscc  ity  on 
the  experimental  results. 

6. 3. 2. 2.  A Posteriori  Evaluation  of  the  Macroscale 

If  it  is  assumed  that  the  maximum  value  of  At^  is  reliably 
inferred  by  the  procedure  described  in  6. 3. 2.1  and  if  it  is 
assumed  that  At^  for  the  short  and  long  time  analyses  are  the 
same,  then  A*  may  be  calculated  from  the  observed  slope  of 
the  experimental  data  for  large  UQt/M  values.  These  calcula- 
tions have  been  executed;  it  is  found  that  the  dimensionless 
slope  (see  E.4.2.)  is  relatively  large:  i.e.,  0.694,  and 
that  the  corresponding  A*  is  also  relatively  large.  Figure  E.3. 
provides  a numerical  value  of  3.9  for  A*;  however,  the  figure 
also  shows  that  a slight  change  in  the  evaluated  slope  will 
create  a large  change  in  the  inferred  A*  value.  Hence,  the 

indicated  A*  value  is  rather  uncertain.  Figure  34  shows  that 
*2  . 

z is  not  particularly  sensitive  to  the  value  of  A*. 


(Figure  E.3.  is  specifically  for  At/At  = 3.5.  However,  the 
magnitude  of  the  slope  will  be  only  slightly  affected  by  this 
ratio;  see  Figure  34. 


7.  SUMMARY  AND  CONCLUSIONS 

A viable  experimental  technique  to  measure  the  instantandous 
concentration  (D  of  a non-dif fusivc  scalar  contaminant  along 
a line  segment  parallel  to  the  principal  concentration  gradient 
has  been  developed.  This  technique,  which  is  discussed  in 
Section  2,  includes  the  capability  to  generate  small  (0.5  urn 
nominal  diameter)  marking  particles  which  serve  as  the  non- 
diffusive  contaminant.  Theoretical  considerations,  developed  in 
Section  4,  show  that  the  particle's  response  time  (to  an  im- 
pulsively applied  motion  of  the  continuum  fluid  element  which 
serves  as  the  host  fluid)  is  sufficiently  small  that  they 
follow  the  continuum  motion  with  satisfactory  accuracy,  see 
Figure  17.  It  is  also  shown  that  the  particle  diffusion  caused 
by  the  Brownian  motion  of  the  host  gas,  has  a negligible  influence 
on  the  inferred  dispersion  properties  of  the  turbulence  field. 

The  experimental  facility  was  a closed  mixing  chamber  fitted 
with  a (horizontal)  splitter  plate  to  separate  the  two  host 
gas  regions  prior  to  their  mixing.  A bi-plane,  square  mesh, 
round  rod  grid  was  used  to  create  the  mixing  field.  Light 
scattering  measurements  from  small  cylindrical  volumes 
(?  ~ 0.36  mm,  d ~ 0.1  mm)  were  used  to  sense  the  presence  or  absenct 
of  the  marking  contaminant.  The  concentration  at  a given 
location  (z)  and  at  a given  time  from  the  passage  of  the  grid 
(t)  will  be  0 or  1 depending  upon  the  origin  of  the  fluid 
element  in  the  scattering  volume.  O'  = 1 if  the  origin  was  in 
the  lower  half-volume.)  Ensemble  averages  of  the  F(z,t)  data 
have  been  constructed  (<T(z,t)>).  The  separate  assumptions  of 
(i)  a constant  turbulent  diffusivity  (D{ ) and  (it)  a Gaussian 
distributuion  of  the  particle  displacement  lead  to  the  same 
analytical  result  that  <r(z,t)>  is  a Gaussian  Cumulative 
Distribution  Function  (C.D.F.).  The  experimental  data  were  fit 
with  such  a C.D.F.  and  the  standard  deviation  of  the  fit  was 
evaluated. 

Conclusion:  The  Gaussian  C.D.F.  provides  a satisfactory 

description  of  the  experimentally  determined 
<T (z,t) The  standard  deviation  of  this 
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distribution  (o  ) provides  an  appropriate  width 

measure  of  the  ensemble  averaged  mixed  region. 

The  experimentally  determined  o and  the  above  noted 

m 

analytical  considerations  allow  to  be  evaluated. 

Conclusion:  The  ratio  of  to  the  host  gas  molecular 

diffusivity  is,  for  air  as  the  host  gas, 
Dfc/D  =73  5 < Uot/M  < 20 

Dfc/D  = 146  20  < UQt/M  < 200 


The  dispersive  property  of  the  decaying  turbulence  field 

. 2 

is  represented  by  the  mean  square  particle  displacement  z ; see 

2 2 P 

(3.2.8.).  [zr  (t)  is  shown  to  be  equal  to  o (t)  for  the  assumed 
p m 

Gaussian  C.D.F.  of  the  particle  displacements).  An  extensive 

2 

theoretical  development  for  the  relationship  of  z^  to  the  Lagrangian 
micro  and  macroscales  which  are  used  to  characterize  the  Lagrangian 
auto  correlation  function,  R , has  been  carried  out  in  Appendix 
D.  This  analysis  has  been  used  to  reevaluate  the  experimental 
results  of  Shlien  [1972],  see  Appendix  E.  The  present  theoretical 
work  utilizes  the  concept  of  a stationary  R (x),  Batchelor  and 
Townsend  [1956],  where  x is  a stretched  time  coordinate 
(x  = C ln[(t  - tQ)/(tr  - tQ)]);  this  feature  is  in  common  with 
prior  analyses.  The  present  theoretical  structure  is  distinct 
in  its  use  of  separate  micro  and  macro  scales  for  the  Lagrangian 
correlation  function  and  for  its  relatively  direct  procedure  to 
relate  the  experimental  and  theoretical  results.  Specific 
conclusions  based  upon  the  comparison  of  the  theoretical  and  ex- 
perimental results  follow. 

Conclusion:  The  experimental  data  of  Shlien  [1972]  are 

well  fit  by  the  present  analytical  formulation. 

The  Lagrangian,  normalized  n icroscale  x*  was 
. 2 

determined  from  the  z data  for  small  values 

P 

of  At/At^;  its  numerical  value  is  0.185.  The 

Lagrangian  normalized  macro  scale  was 

2 

determined  from  the  z^  data  for  large  At/At^ 
values;  its  numerical  value  is  0.213. 


*: 
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Precise  Atr  values  (42  and  98)  are  known  for  the  Shlien  data 
since  the  thermal  marker  was  introduced  with  a heated  wire;  a 
precise  value  of  At  cannot  be  established  for  the  present 
experiment  since  the  mixing  begins  in  the  inhomogeneous  region 
immediately  following  the  passage  of  the  grid.  However,  for 
the  present  study,  the  x*  and  At  evaluations  may  be  jointly 
determined  by  forcing  agreement  between  the  experimental  and 
theoretical  slopes  of  the  z^  vs  At/Atr  curves.  This  procedure 

results  in  a range  of  t*  and  At _ values;  this  range  can  be 

• . r 2 
limited  by  the  following  observation.  z^  (t)  initially  grows 

in  apparent  agreement  with  the  theoretical  description  provided 

by  the  microscale  dominated  region  of  the  R (r)  distribution, 

L 

it  passes  through  a "plateau"  region  where  its  growth  is 
relatively  slow,  and  it  continues  its  arowth  in  apparent  agreement 
with  the  macroscale  dominated  region  of  the  RT  (r)  distribution. 
(Both  the  air  and  Freon  12  host  gases  exhibited  this  behavior.) 

The  initiation  of  the  plateau  allows  a maximum  t*  to  be  defined. 


Conclusion : 


The  maximum  value  of  i*)  is  0.35  and  the 

air 

corresponding  UoAtr/M  value  is  5.17,  see 


section  6 . 3 . 2 . 1 . 

Conclusion:  The  macroscale,  A*,  calculated  from  the 

comparison  of  the  data  to  the  long  time 

theoretical  result,  is  1.9,  from  Section 

6. 3. 2. 7.  This  value  is  recognized  to  be 

quite  uncertain  because  the  calculation 

procedure  shows  that  A*  is  quite  sensitive 

to  small  changes  in  the  At  inferred  from 

the  macroscale  evaluation.  The  "maximum 

value"  of  At  was  used  to  determine  A*. 
r 2 

The  Reynolds  number  dependence  of  the  o and  z~  results  was 
not  expected;  it  was  expected  (d  priori)  that  these  measures 
of  the  dispersion  field  would  solely  reflect  the  large  scale 
motions  and  that  the  scale  of  these  motions  would  be  fixed  by 
the  grid  mesh  size,  M.  However,  the  theoretical  considerations 
of  the  disperion  analysis  provide  a rational  explanation  for  the 
experimental  observations.  Specifically,  it  is  shown  that  i*  may 


be  expected  to  decrease  as  the  Reynolds  number  is  increased  since 

the  constant  C,  in  the  t = C In  (At/Atr)  transformation,  is  Reynolds 

number  dependent.  Hence  the  6 fold  change  in  v between  Freon  12 

2 

and  air  can  be  expected  to  change  z^,  the  mean  square  particle 
displacement. 

Conclusion:  Using  the  air  results  for  the  maximum  t* 

and  the  theoretical  structure  of  Appendix 
D,  the  smaller  dispersion  for  the  Freon  12 
shown  in  Figure  31,  is  made  rational  both 
qualitatively  and  quantitatively  from 
Equations  (6.3.1.)  and  (6.3.3.)  as 


r*  Freon  12,  _ „ t*  Freon  12,  „ „ „ 

— r*Ali } °’69,  ~T*  Air } = 0,74 

experimental  theoretical 

The  availability  of  "instantaneous"  scans  for  F (z,t)  allows 
additional  measures  of  the  mixing  field  to  be  evaluated  as  discussed 
in  Section  5.4.  Statistical  measures  of  the  instantaneous  scans 
are  examined  to  determine  if  the  ensemble  average  mixinq  results 
can  be  traced  to  (1)  a well  mixed  region  for  each  realization  which 
has  the  same  nominal  width  as  the  ensemble  average  field  or  (2) 
a relatively  narrow  mixed  region  which  is  convected  by  the  large 
scale  motions  to  yield  the  observed  width  of  the  ensemble  average 
field.  The  instantaneous  scans  were  shifted  (Az)  to  align  their 
instantaneous  centers.  The  width  reflects  the  magnitude  of 

these  shifts;  reflects  a measure  of  the  Gaussian  C.D.F.  which 
was  fitted  to  the  shifted  and  ensemble  averaged  concentration 
distributions . 

Conclusions:  The  approximate  equality  of  Oj.  and 

(aT  > .8a  ) and  the  small  values  of  o. 

(a.  < .5a  ) imply  that  the  observed  mixinq 

is  strongly  influenced  by  the  relatively 
small  motions  of  the  turbulence  field,  see 
Figure  30. 


Conclusions : 


The  essential  equivalence  of  the  a.  values  for 

AZ 

the  Freon  12  and  air  host  gases  implies  that 
°Az  a measure  of  the  large  scale  motions 

and  that  these  large  scale  motions  are  in- 
sensitive to  the  Reynolds  number. 
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Figure  2.  Detail  of  Turbulence  Generating  Grid 

Note:  Front  wall  and  splitter  plate  have  been  removed  for  this 

photograph 


Figure  3.  Schematic  Representation  of  the 
Acceleration/Deceleration  System 

Notes:  A shaft  connects  to  the  grid-plate  unit 

B hydraulic  oil 

The  steady  state  motion  of  the  piston  represents  balance 
between  the  driving  effects  of  the  compressed  air  and 
the  resisting  effects  of  the  high  pressure  hydraulic 
oil.  The  deceleration  of  the  piston  is  caused  by  the 
pressure  increase  in  the  hydraulic  oil  as  the  conical 
plug  (a)  seats  into  the  tapered  nozzle  (b) ; see 
Appendix  B. 


Fiaure  4.  Velocity  of  Grid  as  a Function  of  Time 


Arqon  Ion 
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Schematic  Representation  of  the  Optical  and  Data  Acquisition  Systems 
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Epmt 


a.  Instantaneous  Scan 

Note:  E.  = threshold  voltage  level  for  the 

comparator  circuit,  see  Figure  13. 


Note:  No  signal  dropouts 


i 

0.2  volts 

T 


Figure  12.  Ink  Tracings  of  Oscilloscope  Photographs  Showing: 

a. )  An  instantaneous  scan  of  the  photomultiplier 

tube  voltage  (E  ) for  the  second  scan 
following  the  pSSsage  of  the  grid 

b. )  A 3-second  time  exposure  of  the  scope  traces 

for  the  fully  mixed  condition  of  large  time 


Figure  14.  Circuit  Diagram  for  the  Comparator  Circuit 
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Figure  15.  Representation  of  the  non-dif fusive  mixinq  field 
Note:  Gas  A = Air,  Gas  li  = Air 

Note:  $ - presence  of  gas  B in  the  scattering  volume;  100  of 

the  205,  data  points  are  shown.  Nominally  every  other 
data  point  in  a scan  is  shown:  ( 1 , 2 , 3 . . . 1 00)  Correspond 
to  2,4,6, . . .200) 

For  the  graph  horizontal  locations:  1 is  the  z value:  +37.1  mm. 

100  is  the  z value:  -35.5  mm 
For  the  graph  vertical  locations:  l * 120  are  the  scan  numbers- 
each  scan  is  executed  in  3.11  msec  and  the  scans  are  repeated 
every  3.6  msec. 

Uc  (3.5  msec)  = 0.865  M. 
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Definition  sketches  for  x (microscale),  x.  (juncture  point)  and  ' 
(macroscale!  for  the  Lagrangian  auto  correlation  function  which  is 
assumed  to  be  stationary  in  the  stretched  x coordinates 
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the  Lagrangian  dispersion  analysis. 


in 


Figure  16b. 


\ \ \\\ti  | >'  no  nil . 

m t of 


«l.  I .» r f r*  view  of  the  ••lxi'u’.  i’u  n,v  mi 

r.\  a| 

, •'<*.1  times  fiv.n  the  * lf;e  of 

the 

grid  .ire 

0 0.02  sec 

il)  O.Oh  !*C 

»ii'  0 . sec 

\v)  1.V6  sec 

v>  ,rC 

.’he  .1 1\  k 

ier  , te  Vs  .It  the  hot  tom  of  \ 

, h. 

esentative  <r(z,t)>  Distributions 
the  Gas  Mixing  Experiments  Presente' 
tandard  and  Probability  Coordinates 


Freon  12 
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apparent  nixed  region  migrated  beyond  the  optically  sampled 
region  of  the  mixing  chamber. 


o 


'J  7 


V 1 

to 

f- 4 

4) 

4) 

• 

E 

H 

• 

3 

S3 

p 

• 

to  H 

E 

g 

0 

04 

p O 

4> 

4> 

4) 

p 

pH 

> 

to 

E 

to 

JQ 

c 

■ <~1 

Tl 

c 

4)  >p 

p 

E 

41 

p 

a> 

0 

E O 

c 

<U 

a» 

(0 

a» 

3 

41 

(0 

4) 

a 

d 

p 

• — i p 

E 

c 

•C 

X 

tp 

O 41 

■H 

4) 

P 

a> 

10 

> xj 

P 

•H 

. — , 

• s 

E 

0> 

E 

E 

~G 

\ 

ON 

4)  3 

G. 

. 0 

0 

4-> 

£ 

H \£* 

-p  C 

X 

p 

p 

•rn 

O 

< 

J A | 

a. 

0) 

4-1 

tp 

•H 

N 

•m 

E 4» 

*» 

x: 

V 

V 1 

id  x: 

x: 

P 

P 

W 

> 

1 

p 

to  P 

4~> 

4> 

4) 

P 

*1 

« 

•m 

i P P 

P 

• 

* 

tp  x> 

C C 

C 

-P 

Cl, 

• 

o c 

4)  4) 

4) 

C 

• 

( 

j p 

ID 

G 

U E 

CJ 

O 

a 

N 

•p 

p 

<T3 

■ H 

o 

• 

'»— * 

<D 

4)  CO  x1 

O 

to  p 

IQ 

u 

O V£>  ft) 

CO 

3 4> 

3 

10 

4-4 

P 

E p 

0 G. 

0 

c 

0 

4) 

3 to  U 

to 

4)  X 

4) 

0 

ID 

J 

C -P  4» 

1 

C 41 

C 

0) 

•H 

G 

E 

i — 1 

0 

fO 

•D 

G 

CO 

0 

3 

h B 41 

to 

0) 

p x: 

P 

(0 

•H 

C 

ID  to 

41 

c 

c P 

C 

-P 

1 

P 

P X) 

P 

<TJ 

id  -t—i 

(0 

c 

<TJ 

P 

O -p  >> 

3 

>p 

p 

P 

<TJ 

o 

• H 

C 

P 3 iP 

to 

c 

to  * 

to 

P 

> 

4> 

r— 4 ■ H 

3 

<T3 

C 4) 

c 

10 

C 

4) 

E 

4-4  p 

4> 

4~> 

-P  C 

•H 

c 

• H 

T3 

■«H 

x: 

E 

to 

•H 

•H 

p 

4-»  x:  p 

c 

4^  -H 

P 

X) 

4) 

-P  4-» 

C 

•H 

0 p 

o 

4> 

o 

P 

a, 

v *H  -*H 

ID 

V 

-C 

p 

— «K 

ID 

X 

C 5 XJ 

U 

-C 

p P 

p 

P 

P 

X) 

4> 

O p 

01 

4-) 

c c 

c 

E 

* 

G 

•P  to  ID 

•H 

41  4» 

0)  4-4 

<TJ 

A 

ID 

II 

P 4) 

to 

E t> 

E 

0 

P 

•H 

P 

iD  E w 

3 

4> 

4> 

o 

to 

•r-i 

P 3 -P 

0 

x: 

t>  4> 

U 

4> 

■ 

N 

C *P 

4> 

4-> 

ID  f' 

D 

rr> 

\ 

* 

4)  O \ 

c 

— i <D 

» — 4 

■D 

| 

N ^ 

to  > -P 

ID 

w. 

P 

P 

•P  M 

V 

o 

4'  O 

P 

0 

to  c 

10 

4> 

o 

O 

p Ip  N 

C 

•p  > 

> 

• 

u 

• 

CN 

O-  o L 

<D 

P 

X)  ID 

T? 

<D 

• 

V 

• 

N 1 

4) 

P 

«D 

• 

• 

P P 

to 

■P 

to 

p 

P 

V| 

41  • 

G 

c 

A 

to 

G 

P 

c 1 

U Jin 

M 

0) 

A 

• H 

4) 

<D 

' 

■p  E <"• 

o 

•H 

..V 

C 

10 

P 3 

•p 

N 

cm 

X) 

0) 

to 

p 

iD  C W 

0 

V 

V 

4) 

P 

G 

4» 

E p 

X 

a 

<TJ 

J 

4>  4> 

to 

1 

1 

•H 

a> 

P 

E 

x:  x:  < 

c 

sc 

p 

H 

g 

O P 

0 

*»—>  *r** 

♦r^» 

G 

03  x> 

•H 

— * 

* 'H 

•H  *H 

4 

G 

CO  o o 

•H 

o u 

u 

(D 

■H  O f— 4 

G 

CM 

CM 

N 

U 

x:  rH  4- 

'<4 

• • 

to 

4j 

to 

ii  x: 

O 

4) 

II 

p p> 

a 

•• 

P 

0 t-pi 

10 

3 

*P 

tn  2 J 

G 

0 

•H 

to 

ID 

4) 

ON 

CM 

4-' 

SC 

•H 

•• 

0) 

C 

si 

to 

P 

•P 

p 

4) 

P 

ip 

X) 

P 

vT* 

0) 

■p 

O 

*H 

Q 

s 

CP 

Figure  30b.  Measures  of  the  Mixing  Region  width  as  determined  from  the  conventional  and  the 
Instantaneous  Scans ...  Freon  12  Carrier  Gas 
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element  of  its  origin,  it  is  inferred  that  It  faithfully  marks  the  lower  half- 
volume fluid  over  the  duration  of  the  experiment  (0.6  sec).  The  7.4  cm  scan 
length  covered  2.92  mesh  lengths  of  the  biplane,  round  rod,  turbulence  gener- 
ating arid.  Air  was  used  for  69  experiments;  Freon  12  was  used  for  100  experi- 
rents  [»„r  ■ 5.8vFre0„  n). 

A principal  result  of  this  investigation  is  that  the  particle  dispersion  is 
dependent  upon  the  mesh  Reynolds  number  (U0M/v);  the  width  of  the  mixing  field 
was  found  to  decrease  as  the  Reynolds  number  was  increased.  A theoretical 
model  of  the  dispersion  was  developed  using  the  suggestion  by  Batchelor  and 
Townsend  [1956]  that  the  Lagrangian  velocity  auto  correlation  function  can  be 
made  stationary,  in  a decaying  turbulence  field,  by  an  appropriate  stretching 
of  the  time  coordinate.  The  present  model  utilizes  both  the  Langrangian  micro 
and  macro  scales  to  describe  this  correlation  function;  numberical  values  for 
these  two  scales  are  inferred  from  the  experimental  data.  The  observed  de- 
crease in  the  particle  dispersion  with  increasing  Reynolds  number  is  supported 
by  a previously  suggested  relationship  between  the  Lagrangian  micro  scale  and 
the  mesh  Reynolds  number,  a new  value  for  the  coefficient  in  this  relationship 
is  proposed. 

Analytical  considerations  are  presented  to  support  the  use  of  a (spatial) 
Gaussian  cumulative  distribution  function  for  the( description  of  the  ensemble 

average  concentration;  viz,  r(z,t)  * [2n  om(t)2]"'i  / exp  - [(a  - <z>)/ 

- o 

2o|n(t)‘:]da.  Satisfactory  agreement  between  this  distribution  and  the  experi- 
mental data  was  obtained;  nm(t)2  was  used  to  evaluate  the  mean  square  particle 
dispersion,  Zp2 . Two  properties  of  the  instantaneous  scans  were  defined:  i)  an 
ensemble  average  measure  of  the  instantaneous  width  of  the  mixing  field,  Oj, 
and  an  ensemble  average  measure  of  the  migration  of  the  instantaneous 
center  lines,  oa2.  The  Instantaneous  scans  are  shown  to  have  essentially  the 
same  width  as  that  of  the  ensemble  average  field  a,  > 0.  8om  whereas  o ex- 
hibits very  little  growth  with  respect  to  time  following  the  establishment  of 
its  nominal  value  at  U t/M  = 8. 
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